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simulations of actual charging events observed on Day 87 and Day 89 (1979) of the 
emission, as well as simulations using test environments, were carried out. 

Numerical models for the simulation of particle emitters and detectors, built into 
NASCAP, were used to analyze the operation of these devices on board SCATHA. The 
effect of highly charged surface regions (so-called "hot spots") on the charging 
of neighboring surfaces was examined. A model for effective conductivity within 
a photosheath was incorporated into NASCAP and used to interpret results from the 
on-board electric field experiment (SC-10) . 

Shadowing calculations were carried out for the satellite and a table of effective 
illuminated areas compiled. 


An analytical treatment of the charging of a large object in polar earth orbit was 
carried out in order to obtain a preliminary indication of the response of the 
shuttle orbiter to such an environment. 
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SUMMARY 


This report describes theoretical investigations per- 
formed to simulate the charging response of the SCATHA 
(Satellite Charging AT High Altitude) or P78-2 satellite in 
geosynchronous orbit. The focus for the investigations was 
a detailed three-dimensional model of the SCATHA spacecraft 
which was used with the NASCAP ( NASA Charging Analyzer Pro- 
gram) computer code . 

Two NASCAP models of SCATHA are described, and their 
material properties are characterized in detail. The charging 
response of individual materials and the full SCATHA models 
is described for a series of test environments and actual 
plasma environments observed by SCATHA in earth orbit. The 
results of NASCAP simulations of actual charging events are 
reported. 

The operation of active control and particle detectors 
were also simulated using NASCAP. Analytical and computational 
models for the operation of the electron and ion guns are pre- 
sented. 

An analytical model of the influence of a small region 
of high potential was developed to simulate the effect of "hot 
spots" on a spacecraft, and the results were compared to the 
NASCAP model. Two code enhancements were implemented to im- 
prove the simulation of photosheath effects: one incorporates 

a model of effective surface conductivity in the photosheath, 
and the other allows the self-consistent calculation of the 
space fields in the photosheath for fixed spacecraft poten- 
tials. The latter option was used to aid in the interpreta- 
tion of the SCATHA electric field experiment. 

A table of effective illuminated areas of the SCATHA 
body was prepared as a supplement to the shadowing tables pre- 
viously generated for selected SCATHA experiments. 


An analytical treatment of the charging of a large 
object in polar earth orbit is presented and its implications 
for the operation of the shuttle orbiter are discussed. 
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1. INTRODUCTION 


This report describes a portion of the work performed 
by Systems , Science and Software on Contract NAS3-21762, 
"Additional Application of the NASCAP Code” . The report 
covers work in which the NASCAP computer code was used in 
conjunction with supplementary analytical models to analyze 
the charging effects of the natural space environment on the 
SCATHA spacecraft and to analyze the combined effects of this 
environment and of the charged condition of the spacecraft on 
the scientific instruments of SCATHA. This work is part of a 
continuing series of analyses designed to assist in the in- 
terpretation of the data collected by the SCATHA spacecraft, 
and to validate and verify the NASCAP computer code as a 
modeling tool for analysis of spacecraft charging. The de- 
velopment of a validated code for the analysis of spacecraft 
charging is one of the goals of the joint NASA/Air Force 
Spacecraft Charging Investigation program. 

Much of the material contained in this report was 
originally prepared for monthly progress reports during the 
contract year. This document consolidates those reports and 
includes additional material to provide a unified and compre- 
hensive description of the SCATHA modeling effort. It is 
assumed that the reader is familiar with the capabilities and 
general features of the NASCAP computer code, which is de- 
scribed in detail in References 1-3. Reference 4 provides a 
summary of the capabilities of the NASCAP program. The first 
description of the SCATHA model was presented at the 1978 
Spacecraft Charging Technology Conference , ^ and a summary 
of some of the work discussed herein was presented at the 
1979 Fall AGU Meeting. 

The NASCAP SCATHA model includes specification of the 
detailed geometrical, electrical, and material parameters of 
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the spacecraft* Complete descriptions of two models are pre- 
sented in Chapter 2, along with a summary of the SCATHA sur- 
face material properties. A discussion of the representation 
of the plasma environment used by NASCAP is given in Chapter 3. 

Chapter 4 reports on a number of aspects of surface 
charging. These include the effect of the environment repre- 
sentation, material properties, the angular distribution of 
the incident particle flux and the effect of high energy 
particles on the charging and discharging of dielectrics. 

In Chapters 5 and 6 the charging response of the entire 
SCATHA vehicle is discussed. Chapter 5 reports on the response 
of the four-grid model to a set of fictitious test environ- 
ments, while in Chapter 6 the results of detailed simulations 
of actual charging events, using the one-grid model, on 
Days 187 and 89, 1979, are presented. 

An important feature of the SCATHA vehicle operation is 
the ability to use charged particle beams to control the satel- 
lite potential. In discussions with AFGL experimenters, 
various features of the gun operations were chosen for de- 
tailed study. The results of these investigations are de- 
scribed in Chapter 7. The NASCAP code includes a DETECTOR 
mode which simulates fluxes to specified spacecraft surface 
locations by explicitly generating reverse particle trajec- 
tories. Chapter 8 summarizes the simulation of the response 
of the satellite particle detector experiments using the 
DETECTOR mode. 

A simple analytical model was developed to illustrate 
the range of influence of an isolated highly charging in- 
sulating spot on the nearby spacecraft surfaces, both with 
and without photosheath effects. A description of these 
models is presented in Chapter 9. Photosheath effects away 
from the surface of the vehicle have also been modeled using 
a self-consistent treatment. This work is described in 
Chapter 9. 
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Shadowing tables for selected SCATHA experiments were 
described previously. ^ For analysis of the overall response 
of the vehicle, shadowing of body elements by booms and pro- 
truding experiments must be considered. Chapter 10 describes 
such body shadowing effects. 

A physical model for the charging of objects large com- 
pared to the Debye length of the ambient plasma is developed 
in Chapter 12. The implications of this for the shuttle 
orbiter in polar earth orbit are discussed. 

Finally, the major conclusions of the study are pre- 
sented in Chapter 13, along with a number of recommendations 
for future investigations. 
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2. SCATHA MODEL DESCRIPTION 


The NASCAP program allows the specification of the 
geometrical, material, and electrical properties of the SCATHA 
spacecraft in considerable detail. This chapter describes the 
NASCAP SCATHA models which were used to perform the charging 
calculations described in Chapters 4 and 5. The models are 
similar to a preliminary version described elsewhere. ^ 

2.1 "FOUR-GRID" MODEL 

Perspective views of the four-grid model are shown in 
Figures 2.1 and 2.2. The main body of the satellite is repre- 
sented as a right octagonal cylinder, with the aft cavity 
visible in Figure 2.2. The 11.5 cm grid resolution allows 
the model to reproduce actual SCATHA geometrical features ex- 
tremely well, as shown in Table 2.1. Note in particular that 
the treatment of booms in NASCAP allows the actual boom radii 
to be reproduced exactly. Figure 2.3 illustrates the com- 
putational space in which NASCAP solves Poisson's equation 
for this model. Monopole boundary conditions are imposed on 
the edges of the outermost grid, which is a rectangular prism 
of dimensions 12.8 * 12.8 x 25.6 m. The zone size increases 
by a factor of 2 in each of the four successive grids. This 
doubling of zone size, plus the requirements that booms 
parallel coordinate axes and intercept mesh points in all 
grids effectively force any long booms to pass through the 
center of the innermost mesh. Therefore, the model includes 
only the SC6, SC11, and two SC2 booms, with the orientations 
fixed at right angles to one another. 
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Figure 2.1. Four-grid SCATHA model: side view. The 50 m 

antenna and the SC1-4 boom are not included in 
this model. 



Figure 2.2. Four-grid SCATHA model? bottom view with aft 
cavity visible. 



frtt. 


Computational space surrounding the four-grid 
SCATHA model, showing the nesting of the grids. 
The tic marks along the axes indicate the outer 
grid zone size; the zone size decreases by a 
factor of two in successive grids. 



TABLE 2.1. COMPARISON OF ACTUAL SCATHA GEOMETRICAL 
FEATURES TO FOUR-GRID NASCAP MODEL 

Zone Size - 4.54 in. (11.5 cm) 


Radius 

Height 

Solar Array Height 
Bellyband Height 
SC9-1 Experiment 
SC6-1 Boom 

Surface Area 

Solar Array Area 

Forward Surface 
Area 


SCATHA 

33.6 inches 

68.7 
29 

12.0 

9.2 x 6 x 8 

1.7 (radius) 

118 (length) 

4 

2.16 x 10 sq.in. 
1.23 x 10 4 

0.36 x io 4 


MODEL 

32.0 inches 

68.0 

27.2 
13.6 

9.1 x 4.5 x 9.1 

1.7 

113.2 

4 

2.11 x io sq.in. 
1.15 x 10 4 

0.34 x io 4 
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2.2 


"ONE-GRID" MODEL 


In addition to the four-grid model it is possible to 
represent the SCATHA satellite entirely within grid one, in- 
creasing the zone size to 19.6 cm. All four materials on the 
SSPM's are resolved but the booms are not to scale. However, 
this so-called "one-grid" model shows a similar charging 
response to the more detailed model above but uses considerably 
less computer time. 

For this reason the one-grid model was used for the 
charging simulation using actual data measured in space by 
SCATHA (Chapter 5) . 

NASCAP generated material plots and perspective views 
are shown in Figures 2.5 through 2.13. 

2 . 3 SURFACE MATERIALS 

The models include the specification of 15 distinct 
exposed surface materials, each of which is specified by the 
values of 14 parameters. The surface materials are described 
in Table 2.2. Wherever possible, experimentally measured 
values for all parameters were used; where this has not been 
possible, suitable estimates based on the properties of 
similar materials were used. Table 2.3 gives the values of 
the material parameters which were used in the calculations 
reported herein. 

The exposed materials in the four-grid model are il- 
lustrated in Figure 2.4, in which the locations of several 
SCATHA experiments are also shown. Experiments at the ends 
of the SCATHA booms were modeled as single boom segments whose 
radii were chosen to reproduce the exposed surface area of 
the actual experiment. NASCAP generated plots of exposed 
surface materials are included in Appendix A. 
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Figure 2,4c. Four-grid SCATHA model with exposed surface 
materials illustrated. 
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Figure 2.6. One-grid SCATHA model material plots 
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TABLE 2.2 


EXPOSED SURFACE MATERIALS 


GOLD: 

SOLAR: 

WHITEN: 

SCREEN: 

YELOWC : 
GOLDPD : 

BLACKC : 
KAPTON : 
ASTROQ : 
TEFLON : 
INDOX: 
YGOLDC : 

ML12 : 


ALUM: 

BOOMAT: 


gold plate 

solar cells, coated fused silica 

non-conducting white paint (STM K792) 

SC5 screen material, a conducting 
fictitious material which absorbs but 
does not emit charged particles 

conducting yellow paint 

88 percent gold plate with 12 percent 
conductive black paint (STM K748) in a 
polka dot pattern 

conductive black paint (STM K748) 

kapton 

Si0 2 fabric 

teflon 

indium oxide 

conducting yellow paint (50 percent) 
gold (50 percent) 

ML12-3 and ML12-4 surface, a fictitious 
material whose properties are an average 
of the properties of the several materials 
on the ML12 surfaces 

aluminum plate 

platinum banded kapton 


TABLE 2.3. MATERIAL PROPERTIES FOR EXPOSED SURFACES 
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TABLE 2.3. (Continued) 


a The materials are described in Table 2.2. 
b The thirteen properties are as follows: 


Property 

1: 

Relative dielectric constant for 
insulators (dimensionless) . 

Property 

2: 

Thickness of dielectric film or 
vacuum gap (meters) . 

Property 

3: 

Electrical conductivity (mho/m) . 
The value «° indicates a vacuum gap 
over a conducting surface. 

Property 

4: 

Atomic number (dimensionless) . 

Property 

5: 

Maximum secondary electron yield 
for electron impact at normal inci- 
dence (dimensionless) . 

Property 

6: 

Primary electron energy to produce 
maximum yield at normal incidence 
(keV) . 

Properties 7-10: 

Range for incident electrons. 
Either : 


8 10 
Range * P,E + P^E 

where the range is in angstroms and 
for the energy in keV, 

or 

P 7 * -1. to indicate uss of an 
empirical range formula 

Pg * density (g/cm^) 

P^q » mean atomic weight (dimension- 
less) . 


Property 

11: 

Secondary electron yield for normally 
incident 1 keV protons. 

Property 

12: 

Proton energy to produce maximum 
secondary electron yield (keV) . 

Property 

13: 

Photoelectron yield for normally 
incident sunlight (A/m^) . 

Property 

14: 

Surface resistivity for insulators 
(ohms) . 


c The dielectric constant and thickness for the boom surfaces 
were chosen to reflect the effective capacitance to the under- 
lying cable shield. 
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The model also includes six distinct underlying con- 
ductors: spacecraft ground, the reference band, and the four 

experimental mountings SC2-1, SC2-2, SC6-1, and SC6-2. Each 
of these conductors can be separately biased or floated with 
respect to spacecraft ground, and each conductor is directly 
capacitively coupled to spacecraft ground. The values em- 
ployed for these capacitive couplings are given in Table 2.4: 
these values were chosen to represent the capacitance of di- 
electric spacers separating the conductors from ground. 

To improve agreement between NASCAP predictions and 
experimental observations, the values of the material para- 
meters are continually being updated, as better information 
becomes available. In particular some effort has been ap- 
plied to the problem of secondary emission due to the impact 
of both electrons and ions at the spacecraft surface. 


TABLE 2.4. CAPACITIVE COUPLINGS EMPLOYED WITH 
SCATHA MODEL 

Conductor Capacitance to Ground (pf) 


2: 

SC2-1 

30 

3: 

SC2-2 

30 

4 : 

SC6-1 

240 

5: 

SC6-2 

30 

6 : 

Reference Band 

250 
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2.4 ELECTRON-INDUCED SECONDARY EMISSION 

The formulation of secondary emission due to incoming 
electrons has been improved by the use of Letter stopping 
power data. The secondary electron yield 6 for a primary 
beam normally incident is directly proportional to the 
stopping power S Q at tne incident energy. 

R 

6 = CS Q / f (x) e“ ax dx 

0 

where f(xj <?e; cribes the change in stopping £ wer as a 
function cf depth x, and R is the "range" of the primary 
electron . 

By fitting S Q to ab-initio calculations by Ashley, 
et al. much more physical values of 6 were obtained for 

materials like kapton, especially at energies above 10 keV 
where they were previously too high (Figure 2.14). 

The stopping power is fit in such a way as to force 
the range into a bi-exponential form. 

n l n 2 

R (E) = bjE x + b 2 R 

S q (E) = ( n^E 1 + n 2 b 2 E / 

The parameters n^b^n 2 b 2 are given in Table 2.5 for a number 
of materials. 
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E (Energy, kev) 

Figure 2.14. Comparison of updated and original secondary yields for kapton 








TABLE 2.5. FOUR PARAMETER STOPPING POWER FITS 


s * 

n l’ 1 

n^b^E + n 2 b 2 

n,-l “ 1 
E 


Material 

h 


^2 

^2 

Au 

88.79 

0.92 

53.48 

1.73 

Ag 

84.46 

0.82 

79.43 

1.74 

A1 

153.7 

0.80 

220.0 

1.76 

Si0 2 

116.3 

0.81 

183.1 

1.86 

Polystyrene 

88.43 

0.67 

371.7 

1.79 

Kapton 

71.48 

0.60 

312.1 

1.77 

Teflon 

45.37 

0.40 

217.6 

1.77 

Solar* 

77.50 

0.45 

156.1 

1.73 

♦Stopping power of Al 2 0 3 

used. 



ION- INDUCED 

SECONDARY EMISSION 


The NASCAP model assumes 

that 

all positively charged 


species in the plasma environment are protons H + . Measure- 
ments made in geosynchronous orbit indicate in fact that 
often up to 80 percent of the ions present are 0 + rather than 
H + . This observation has called into question the ion-impact 
induced, secondary electron current, calculated by the code 
assuming a purely proton environment. 

Secondary emission of electrons following ion-surface 
impact can occur via two mechanisms. 

a. Potential Emission 

This occurs via transfer of ion potential energy to 
lattice electrons at metal surfaces. Electrons tunnel into 
the potential well formed by the adsorption of the ion on 
the surface, neutralizing the ion, which then auto-ionizes. 
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It is a low energy phenomenon (<20 eV) and is unimportant in 
most of the energy regime associated with NASCAP (0-50 keV) . 

b. Kinetic Emission 

Here emission results from the direct transfer of ion- 
kinetic energy to lattice electrons, and depends, in a compli- 
cated way, on the projectile/target atomic collision cross- 
section. It is the dominant mechanism in the energy range of 
interest. 

The yield for both mechanisms does not appear to depend 
in any predictable way upon atomic number. The mechanism for 
potential emission is almost chemical in nature and depends 
much more upon electronic structure than nuclear mass. The 
most identifiable trend appears -o be an increasing yield 
for projectile ions having greater electron affinities. 

While the collision cross-section central to kinetic 
emission increases with atomic number of the projectile ion, 
the yield of escaping secondary electrons involves a trade- 
off between factors such as the efficiency of energy transfer 
per collision and the depth of penetration of the ion. This 
is rather poorly understood and experimental studies with 
rare gas ions impinging on clean metal surfaces show an ir- 
regular dependence of yield upon atomic number. For example, 
the energy-yield curves for Kr and Xe incident on Mo cross 
twice within the range 6-10 keV. 

Conclusions 

A table of secondary yields for ions of 1 keV incident 
on Mo is shown below. 

Ion Yield (at 1 keV) 

H + 0.250 

He + 0.276 

0 + 0.178 
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In this case the yields for H + and 0* are of a similar magni- 
tude. Very little additional data is available. The data 
that is available is often subject to large errors because of 
the sensitivity of measurements to the nature of the test 
surface. This coupled with the fairly large uncertainties 
in the measurement of 0 + /H + ratios in space at any particular 
time leads us to the conclusion that adjustment of the code 
and/or data to take the presence of 0 into account is not 
justified. The magnitude of the adjustments to be made are 
smaller than the additional uncertainties that would be 
introduced . 


3. REPRESENTATION OF THE PLASMA ENVIRONMENT 


i 

i 


r 

I 


The NASCAP code recognizes plasma spectra in three 
forms: Maxwellian, double Maxwellian and a tabular form as 

described in Appendix C. The tabular data is integrated 
using the DIRECT mode. 

Measurements made on board SCATHA have provided in- 
formation about the spectrum of the plasma environment in 

the form of tabulations of the distribution function, for 

2 5 

both electrons and ions, in the energy range 10 to 10 eV. 
This data is based on observations of energy flux, made by 
the SC9 detector, averaged over a 16 second period, and so has 
certain associated limitations; i.e., 

a. When the satellite is charged the spacecraft poten- 
tial affects the energy of charged particles reaching the 
surface. For example if the spacecraft potential is -2000 V, 
protons are attracted to it and their energy is increased by 
2000 eV at the surface, compared to infinity. In the same 
way electrons are repelled so that their energy is reduced 
by 2000 eV at the surface, and those electrons with initially 
less than 2000 eV of energy do not reach the surface at all. 

This distorts the distribution function observed at 
the surface. While the shape of the distribution function 
f(E> is unaffected by spacecraft potential, the energies 
associated with each value are affected. The particle re- 
pelled by the potential <+> has q x <j> less energy at the sur- 
face than it did at infinity (q is the charge on the parti- 
cle) . Hence instead of associating the value of the distri- 
bution function f 00 (E) with energy E, the value is associated 
with particles of energy E - q<J> , i.e., 

fJE) = f g (E - q<M 
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Knowing the spacecraft potential, this shift in ener- 
gies is easily corrected for, providing the repelled species 
can reach the surface at all. However, those particles with 
energies at infinity less than q$ cannot and so there is no 
distribution function information measured at these energies. 

The lack of information for the repelled species 
(almost always electrons) , below the spacecraft potential, 
represents a major limitation of the data as measured. This 
limitation does not apply to the attracted species (ions) 
since particles of all energies reach the spacecraft. 

b. The energy flux is weighted by the energy 

<30 

<Energy Flux> J E 2 f(E) dE 
0 

Measurements at low energies have a lower intensity 
and poorer signal to noise ratio. Hence the estimates of 
distribution function derived from these measurements are 
less reliable at the lowest energies. 

NASCAP requires spectral information in the form of 
a distribution function at infinity. This presents no dif- 
ficulty for the attracted species but, for reasons discussed 
above, can lead to complications for the repelled species. 

3 . 1 DIRECT INTEGRATION 

In particular, the absence of spectral information 
concerning the repelled species for energies below the 
spacecraft potential has prevented a successful simulation 
involving direct integration of the observed data points. 

In principle, the missing data could be replaced by an 
extrapolation of the known data. This requires some assump- 
tion about its functional form and a fit to this form. Two 
approaches were tried. 
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a. The whole of the known data was fit to a 
Maxwellian form. 

b. The known data of lowest energy was fit to a 
Maxwellian form. 

The former approach was unsuccessful because the effective 
temperature of a Maxwellian fit increases with energy (i.e., 
the known data has a high energy non -Maxwellian "tail”). 

Thus the fit had a higher temperature than that appropriate 
for the extrapolated energy range. Figure 3.1 shows this as 
Case "C" . The higher temperature leads to a lower density 
than would be expected by smooth extrapolation of the known 
data (Case "B") , and hence to anamolously low values for 
the distribution function in the extrapolated region. Such 
a distribution leads to false multiple equilibrium potentials. 

The second approach, fitting only the data at energies 
close to the spacecraft potential, can lead to unphysically 
high values for the density and hence anamolously high values 
of the distribution function in the extrapolated region 
(Case "A”) . This problem arises because the fit is based on 
an unrepresentatively small sample of points. 

Stable predictions were obtained using direct inte- 
gration of the data when only the high energy points were 
used, and all of the remainder replaced by a fit to these 
points. Under these circumstances, it is just as reasonable 
to replace the whole spectrum with a suitable analytic 
representation instead. 
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3.2 


PITTING OP THE DATA 


Two types of fit to the data were made. 

a. Single Maxwellian — fitted by moments 
(Appendix D) . 

b. Double Maxwellian - fitted by least squares 
analysis. 

In general , the double Maxwellian fits were better able to 
represent the non-Maxwellian character of the data (Figures 
3.2 and 3.3). Even though the fits were good from 100 eV 
to 100,000 eV, the low energy components of the double 
Maxwellians were often unphysical with very high densities 
and low temperatures. This was particularly true for the 
electrons in cases where the satellite was charged to several 
thousand volts negative. Under these conditions only the 
very end of a Maxwellian component with a temperature, say 
one-tenth the spacecraft potential, contributes to the flux 
at the surface. Such contributions can only be considered 
as noise. 

Increasing the "cutoff", below which data points are 
ignored, to 1000 eV greatly improves the physical picture 
provided by the fit. This, coupled with an enforced lowest 
value of at least half the spacecraft potential for the 
temperature, leads to a series of double Maxwellian fits 
that were both physical and accurate in all energy regions 
(Table 3.1) . 
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3.3 CHARGE NEUTRALIZATION 

While the fits accurately represent both the ion and 
electron data points as observed, the electron densities tend 
to be as much as a factor of ten higher than the ion densi- 
ties at the same time. This implies that the plasma surround- 
ing the spacecraft is highly non-neutral. This unphysical 
result may be due to a systematic error in the measurement 
of energy flux for one of the species. Independent measure- 
ments indicate that the ion densities are the best choice 
for the correct value. To neutralize the plasma we renormalize 
the electron densities so that they are equal to the ion den- 
sities. This would be a simple matter if the ions were all 
protons (as the conversion to distribution function assumes) . 
However, measurements by Kaye, et al. show that 0 + is 

often the dominant species. This does not affect the fluxes 
as calculated by NASCAP because the code also assumes all the 
ions are protons and so the error is cancelled out. However, 
a factor of (mass) ' is carried over into the values of the 
distribution function and hence the estimate of density N. 

For a pure 0 + environment an electrically neutral plasma 
— 1/2 + 

would have N * (16) N as calculated assuming that all 
ions were protons; i.e., 

N + = 0.25 N - 

If only a fraction a are oxygen 

N + = [0.25a + (1-a) ] N - 
N~ = N + /(l - 0.75a) . 

To correct the values of the density for the electrons the 
values of N~ and N2 obtained by the fitting procedure were 
multiplied by the factor f. 
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Hence we arrive at the final representation of the environ- 
ment, shown for Day 87, 1979, in Table 3.1. 
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4. SCATHA MATERIALS CHARGING RESPONSE 


To examine the sensitivity of the charging response to 
the method of fitting, and compare the fits with directly inte- 
grated tabulated data, the NASCAP support code, MATCHG, was 
used to study just one material. The material chosen was 
"SOLAR", the solar cell cover glass that forms most of the 
exposed surface area of the SCATHA satellite. The silica 
cover glass is coated with a non-reflective MgFj layer and we 
assume that is has the same material properties (e.g., second- 
ary emission yield, etc.) as MgFj. The MATCHG predicted 
equilibrium potentials for a sphere covered with "SOLAR" 
charging under the influence of the three different repre- 
sentations of the environment are compared in Table 4.1. The 
parameters associated with these representations are summarized 
in Tables 4.2 and 4.3. 

TABLE 4.1. CHARGING RESPONSE (kV) OF SOLAR AS A FUNCTION 
OF ENVIRONMENT REPRESENTATION 

59873 


Environment 

59813 

59873 

61073 

(Neutral) 

Single 

Maxwellian 

-14.5 

-22.3 

- 17.2 

f- 1 
• 

<N 

H 

1 

Double 

Maxwellian 

-17.6 

-20.3 

-18.9 

- 

Direct 

Integration 

-16.2 

-22.6 

-17.9 

- 

Satellite 

Potential 

- 1.9 

- 7.1 

- 4.8 

- 7.1 
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TABLE 4.2. SINGLE MAXWELLIAN FITS TO DAY 87 MEASURED 
DISTRIBUTION FUNCTIONS 


Time 

(sec) 


Species 

Densit 
(10 6 nr 

y Temperature 

3) (keV) 

59813 

59813 

(uncorrected) 

Electrons 

Ions 

0.79 

0.086 


8.7 

12.0 

59873 

59873 

(uncorrected) 

Electrons 

Ions 

0.98 

0.15 


12.0 

9.9 

61073 

61073 

(uncorrected) 

Electrons 

Ions 

0.95 

0.20 


11.0 

12.0 

59873 

59873 

corrected 
for 0+ 

Electrons 

Ions 

0.28 

0.15 


12.0 

9.9 

TABLE 4,3. DOUBLE MAXWELLIAN FITS TO DAY 87 MEASURED 

DISTRIBUTION FUNCTIONS (UNCORRECTED FOR 0 + ) 

Time 

(sec) 

Species 

Density Temperature 

(10 6 m~J) (keV) 

Density 
(1C 6 m" 3) 

Temperature 

(keV) 

59813 

Electrons 

0.29 

5.9 

0.48 

11.0 

59813 

Ions 

0.020 

0.70 

0.043 

27.0 

59873 

Electrons 

0.062 

4.7 

0.87 

12.0 

59873 

Ions 

0.037 

1.0 

0.096 

14.0 

61073 

Electrons 

0.64 

5.1 

0.45 

16.0 

610 '3 

Ions 

0.084 

2.2 

0.11 

24.0 
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Agreement between the two fits and the directly inte- 
grated spectra are very good. The direct integration is only 
possible because a single Maxwellian fit to the known data is 
used to fill in the part cut out by the spacecraft potential. 
Furthermore, a much more stable result is obtained when in 
addition the first 2000 eV of known data for both the ions 
and electrons is replaced by the Maxwellian fit. The points 
replaced are highly non-Maxwellian. This could be a real 
phenomenon or it may be due to the poor signal/noise ratio 
inherent in measurements of energy flux at low energies. In 
either case, including it leads to erratic changes in the 
calculated net current and prediction of more than one equi- 
librium potential. These additional predicted potentials 
occur at unphysically low values and are almost certainly 
false. 


Included in Table 4.1 are the spacecraft ground poten- 
tials observed on the SCATHA satellite when the distribution 
function data was being measured. MATCHG predicts considerably 
more charging than observed. There are at least three reasons 
for this: 

1. The satellite is not a sphere covered in SOLAR. 

2. In space the charging environment is not constant, 
particularly at the times when our test spectra 
were measured. MATCHG assumes a steady charging 
environment . 

3. The measured distribution functions are highly 
non-neutral (see Tables 4.2 and 4.3). This is 
rather unphysical and as explained in Chapter 3, 
indicates a systematic deficiency in the measure- 
ment in the density of at least one of the species. 

The equilibrium potential predicted with the environment cor- 
rected for the presence of G + is much closer to the observed 
satellite potential. 
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4.1 


THRESHOLD EFFECT 


The three envii onments chosen all charge the SOLAR 
material to high potentials in all representations. For 
environments with only a small tendency to charge the material, 
very small changes in the representation of the plasma can 
lead to qualitatively very different charging predictions. 

Figure 4.1 shows the effect of varying the electron 
temperature parameter in the 59873 single Maxwellian repre- 
sentation. There is a definite threshold for charging that 
occurs at T * 6.9 keV. Below this temperature, no charging 

is predicted, while only 0.3 keV above it, at T * 7.2 keV, a 

* 

potential of 4.3 keV is forecast. Such small changes in 
temperature are well within experimental error for measure- 
ment and fitting procedures. For "borderline" environments 
close to the charging threshold quite different predictions 
are possible according to the representation employed. 

A similar phenomenon is presented when we examine the 
effect of changing important material parameters on charging. 
Figure 4.2 shows a threshold for charging in the 59873 en- 
vironment (as given in Table 4.2), when the maximum secondary 

electron yield 6 * 3.1. Above this value the net current 

max 

is positive at zero potential and the spacecraft never 
negatively charges, while below 3.1 rapid charging to several 
keV negative is predicted. The usual value used by MATCHG 
of 2.05 is well into the stable charging range. However, in 
borderline cases small changes and errors in material proper- 
ties can lead to significantly different predicted potentials. 

Another factor that can affect charging is the angular 
distribution of the incident fluxes. 


The particularly high value of 4 . 3 keV predicted here is due 
to the non-neutrality of the plasma. 
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max 

Figure 4.2. Charging characteristics of a solar sphere as a function of 
secondary emission yield 

max 



4.2 


ANISOTROPIC FLUX 


The SCS detector on the "belly band" of the SCATHA 
satellite measures the angular distribution of the energy 
flux of the surrounding plasma in the plane of satellite 
rotation (Figure 4.3). For an isotropic plasma the flux is 
constant over a rotation giving a circular radial plot (Fig- 
ure 4.4). For an anisotropic flux (i.e., one having particles 
aligned preferentially in one direction) the measured plot is 
distorted. The degree of distortion increases with the degree 
of anisotropy. 

The most extreme case consists of a narrow beam lying 
in the plane of rotation. As the satellite rotates, the angle 
of incidence of the beam, 0 Q , oscillates between 0 and ir/2. 

To investigate the effect of anisotropy on spacecraft charging, 
we simulated this situation using MATCHG to calculate the 
potential of an aluminum plate under the influence of a beam 
with an oscillating angle of incidence 0 • 

The secondary emission yields and the backscatter all 
increase with increasing angle of incidence, and we expect 
the aluminum to be driven more positive at high angles and 
more negative at low angles of incidence. Figure 4.5 shows 
a plot of potential against time (angle) confirming our ex- 
pectations . 

When the charging response is fast compared to the 
period of the oscillations, we see the potential of the plate 
oscillate in time with the angle of incidence of the beam. 

The faster the charging response (or slower the rotation) , 
the greater the amplitude of the oscillations. For a very 
fast oscillation, the plate would respond only to an average 
environment and no oscillations would be seen. The initial 
oscillations are erratic and more pronounced before settling 
into a regular pattern. This behavior reflects the two in- 
fluences on the net current to the plate; potential and 0 Q . 
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Figure 4.4. Representation of the angular distribution in the 
rotation plane, r represents the magnitude of 
the flux at angle 6. 


50 







Initially, these may be out of phase leading to transient 
potentials that are too high or low. 

The mean of the oscillating potential is not the same 
as that for an isotropic plasma of the same density and tem- 
perature. This is because of different incident currents 
arising out of distributions with the same flux normalization. 
A directional beam lying in the rotation plane has but one 
angle of incidence, 0 Q , so the incident current is proportional 
to cose. The mean incident current is thus the mean of cos6. 


Beam: 


r B «/ 


Tt/2 


tt/2 


cos0 de/ 


/ 


d e - 


For an isotropic distribution an angle of incidence 0 Q in 
the rotation plane has associated with it all other angles 
of incidence due to particles arriving from above and below 
it. The average of these is sine. Thus the current from an 
isotropic distribution is proportional to cos0 and sin0, and 
the mean is the mean of cos0 sin9: 


it/2 it/2 

Isotropic: ql J* cos6 sin0 d0/ f sin0 d0 ■ j 

0 0 

Thus, the beam current exceeds the isotropic by a factor of 
4/tt x 1 . 27 . 

This illustrates some important points regarding the 
measurement of flux distributions in space: 

1. If a detector measures an average flux over a 
rotation and assumes that it arises from an 
isotropic plasma, the actual current will be 
underestimated by an amount that will increase 
with the increasing directionality of the true 
angular distribution (reaching a maximum of 4 /t 
for a beam) . 
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2 . 


If a detector measures actual average current , 
then in the same way flux (density) will be over - 
estimated. 

3. If a detector measures the angular distribution 

of the flux, information in both the perpendicular 
and parallel directions must be known, or implied 
to infer densities and currents. 

4 . For a "loss-cone" (negative aligned component) 
the reverse of 1 and 2 apply. 

4.3 RADIATION- INDUCED BULK CONDUCTIVITY 

In addition to the electron flux with energy below 
100 keV, electrons with energies up to 5000 keV have been 
observed by detectors on board SCATHA. This high energy 
radiation makes an insignificant contribution to the total 
incident electron current but nevertheless can influence dif- 
ferential charging of insulators on a spacecraft 

When high energy radiation, such as a 300 keV electron, 

passes through an insulator such as kapton, electrons can be 

promoted into the normally empty conduction bands and increase 

r 13 * 

the bulk conductivity o. Frederickson 1 J has represented 
this by the equation 

a - <D + o 

o 

where D is the radiation dose rate and < is the coefficient 
that depends upon the nature of the material. a Q is the con- 
ductivity in the absence of radiation. As the flux and hence 
dose rate increases, the radiation-induced conductivity in- 
creases. For a sufficiently high flux this could limit the 
potential differences that can build up between an insulator 
and the underlying conductor. 

To investigate this question we use MATCHG to predict 

-4 

the potential of 0.005 inches (1.27 * 10 m) thick kapton 
subject to the single Maxwellian representation of the 59873 
environment, with a range of values for the bulk conductivity. 
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The results are shown in Table 4.4. The fluxes corresponding 

to each value of a can be estimated using an experimental 

r 14 1 

result of Treadaway et al. He found that a 0.002 inch 

-4 -2 

(5.08 x 10 m) kapton film subjected to a 0.05 pA cm beam 

of 300 keV electrons accompanied by a 0.2 nA cm beam of 

10 keV electrons charged to -1600 (±300) V. Simulating this 

-14 

experiment with MATCHG implied a value of 4.67 * 10 mhos 
m” 1 for the bulk conductivity a. Assuming that o Q is 
insignificant we can estimate k. 

a = kD 

• —8 —2 
The dose rate D arises from 5 x 10 Am of 300 keV 

4 -2 -1 -1 

electrons, i.e., a flux of 1.67 x io electrons cm s sr 
keV This is equivalent to a dose rate of 1.2 rads s" 1 . 

-14 -1 -1 

k = 4 . 0 10 mhos m rad s 

This value is rather higher than Frederickson ' s estimate of 
10” 15 -10" 16 mhos m" 1 rad" 1 s. f151 

As we can see from Table 4.4, as soon as the dose 
2 -2 -1 -1 

reaches ^10 electrons cm s keV there is a significant 

drop in the potentia] difference that the kapton film can 

support. Since a 0.005 inch la' f er of kapton is typical of 

the insulating materials found on satellites, this result 

2 

suggests that in environments with doses higher than 10 
electrons cm * s sr keV , the radiation induced con- 
ductivity may play a significant role in preventing acute 
differential charging and hence discharges. 

Figure 4.6 shows a plot of the data in Table 4.4. The 
vertical lines are drawn to represent the typical values for 
200 keV flux on days 146, 87 and 114. [161 Days 146 and 114 
are examples of the lowest and highest extremes documented 
so far = We see that fluxes in the range where radiation 
induced conductivity appears to be important are common. 
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It will be interesting to discover, as more data becomes 
available, if there is any correlation between the hig>* 
energy flux and discharges on board SCATHA. 


TABLE 4.4. THE EFFECT OF RADIATION -INDUCED CONDUCTIVITY ON 

THE CHARGING OF 0.005 INCH (1.27 x 10**4 m) KAPTON 
FILM AS PREDICTED BY MATCHG 


300 keV 
Incident 
Current 
-2 

pA cm 

Differential 
Flux (F) 
Electrons cm“2 

s” 1 sr” 1 kev" 1 

Conductivity 

a 

* 

Potential 

Volts 

0 

0 


0 


-15500 

0.003 

1.0 x 

10 1 

2.8 x 

10- 17 

-15300 

0.03 

1.0 x 

10 2 

2.8 x 

10- 16 

-13700 

0.05 

1.67 x 

10 2 

4.67 x 

10- 16 

-12800 

0.5 

1.67 x 

10 3 

4.67 x 

10- 15 

- 5600 

5.0 

1.67 x 

10 4 

4.67 x 

lo’ 14 

- 1000 

50.0 

1.67 x 

10 5 

4.67 x 

!0- 13 

- 100 

500.0 

1.67 x 

10 6 

4.67 x 

xo- 12 

- 0 

* 
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Penotratinci Differential Electron Flux (electrons cm s sr keV ) 

Figure 4.6. Plot of equilibrium potential of 5 mil kapton coating a grounded con 
ductor versus bulk conductivity. 





5. SCATHA CHARGING SIMULATIONS USING TEST DATA 


This chapter comprises descriptions of SCATHA charging 

simulations in several typical, but fictitious, environmental 

* 

condit. ons at geosynchronous orbit. Charging in eclipse 
using the "hiqh" and "moderate" double Maxwellian plasmas is 
described in Section 5.1 and 5.2, respectively. Sunlight 
charging simulations are discussed in Section 5.3. The in- 
duced charging event on day 89 is of special interest due to 
the SC2 failure and telemetry upset which occurred; modeling 
of the SCATHA charging response during this event is discus- 
sed in Section 5.4. 

The four grid SCATHA model as described in Chapter 2 
was used for all of the simulations reported here. The 
spherical probe current collection model and the "NORMAL" 
secondary yield formulation were used throughout. Monopole 
boundary conditions were imposed at the outer (fourth) grid 
boundary, and simulations were begun with the spacecraft un- 
charged, unless otherwise indicated. 

5.1 ECLIPSE CHARGING IN HIGH TEMPERATURE ENVIRONMENT 

The charging simulation described in this section used 
a high temperature double Maxwellian ambient plasma model 
given in Table 5.1. The spacecraft response followed a pattern 
typical for a highly charging environment: rapid overall 

charging to several kilovolts followed by much slower develop- 
ment of differential charging. Figures 5.1 and 5.2 illustrate 
the potential contours after three seconds of charging, when 
the overall charging to -6.4 kV is complete. Differential 
charging was then followed for 3200 seconds: Figures 5.3 

and 5.4 illustrate the final potential contours. 

These computations were carried out early in the program, 
prior to the updating of material properties shown in 
Table 2.3. The conclusions remain valid, however. 
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SCATHA IN ECLIPSE; S/C GROUND = -6200 V; CONTOUR STEPS = 500 V. 

Figure 5.2. SCATHA potential contours. Timt = 3 seconds. *1 
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SCATHA IN ECLIPSE; S/C GROUND = -6400 V; CONTOUR STEPS = 500 V 


Figure 5.4. SCATHA potential contours. Time = 3230 seconds. 





Two types of differential charging occur during this 
period. The first involves development of electric stresses 
internal to insulating materials whose surfaces charge rela- 
tive to underlying conductors. Maximum stresses of this type 
were 3 * 10® V/m for TEFLON insulating patches on the belly- 
band; far below breakdown thresholds for this material. Dif- 
ferential charging along adjacent satellite surfaces can also 
occur. In this simulation, the segments at the ends of the 
SC2 and SC6 booms representing experiment mountings were al- 
lowed to float with respect to spacecraft ground. As a result, 
differential charging of surfaces near these boom ends becomes 
severe, as the SC6 tip (GOLD) reaches -3.5 kV and the SC2 tips 
(BLACKC) charge to -9 kV. (The SC11 boom remains nearly an 
equipotential since the magnetometer experiment at its end is 
not distinguished in the model.) The large surface potential 
gradients near boom tips are clearly visible in Figure 5.4: 
these are the most likely sites for surface flashovers identi- 
fied in this simulation. The development of differential 
charging during the simulation is illustrated graphically in 
Figure 5.5. 

5.2 ECLIPSE CHARGING IN MODERATE TEMPERATURE ENVIRONMENT 

The charging simulation described in this section used 
the moderate temperature double Maxwellian ambient plasma model 
given in Table 5.1. The expected spacecraft response is quite 
different from that described in the preceding section: 
little overall charging is expected since equilibrium poten- 
tials for much of the surface material (GOLD, GOLDPD, SOLAR) 
are near zero. The satellite surfaces charge differentially 
(hence slowly) from the outset: potential contours after 

2300 seconds are shown in Figures 5.6 through 5.8. Even 
though the charging environment is not severe, local charge 
buildup near some surfaces can lead to large differential 
charging, such as near the KAPTON surface of the SCI experi- 
ment in Figure 5.6, and along the SC2 booms in Figure 5.7. 
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Figure 5.5. Potentials versus time for SCATHA model in high 
temperature ambient plasma. 
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Figure 5.6. 


SCATHA potential contours. 
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Figure 5.7. 


SCATHA potential contours. 
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Figure 5.8 


SCATHA potential contours 



Charging of the non-conducting white paint (WHITEN) covering 
the aft surface has important consequences. Fringing fields 
from the aft surface result in the suppression of low energy 
emission from neighboring solar cells and the cavity, so that 
the ground conductor has charged to -143 volts at the time 
shown. The WHITEN surface has reached -390 volts; further 
charging is limited by bulk conductivity. The dipolar charac- 
ter of the potential is clearly displayed in Figure 5.8. The 
reference band conductor (GOLD) reached a potential of -213 
volts, nearer the WHITEN potential than the ground conductor 
due to its proximity to the aft surface. The drift of the 
reference band away from plasma ground due to fringing fields 
will obviously have an impact on SSPM measurements which at- 
tempt to use the reference band as a measure of zero potential; 
this point is discussed further in the next section. 

5.3 SUNLIGHT CHARGING SIMULATIONS 

Two simulations of charging response in sunlight are 
described in this section. Since the differential charging 
timescale for most SCATHA materials is typically 'vlO'* seconds, 
the satellite reaches equilibrium in response to solar illumi- 
nation averaged over many rotations, and the SPINNER model in 
NASCAP is appropriate. In this mode, both a sun direction 
and a spin axis are specified, and the photocurrent expected 
from each surface cell is calculated on the basis of the 
averaged applied solar illumination. Although the relatively 
fast response of instruments such as the SSPM experiments 
could rot be correctly modeled in this fashion, the SPINNER 
mode does serve «rt a convenient and economical procedure for 
studying ex. -„ed equilibrium charging response. 

The results of the first set of calculations are il- 
lustrated in Figures 5.9 and 5.10. These potential contours 
result from charging in the high temperature single 
Maxwellian environment, using the SPINNER model in sunlight; 
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over 5000 seconds have elapsed in the simulation. The sun 
direction was taken to be perpendicular to the SCATHA spin 
axis, so that only the forward, aft, and cavity surfaces re- 
main in darkness. Fringing fields from the aft surface are 
not sufficient to cause charging of the ground conductor: 
it remained at +4 volts throughout the simulation. All the 
boom surfaces also stayed at small positive potentials. The 
charge accumulation of the WHITEN aft surfaces again leads 
to an overall dipolar field, as shown in Figures 5.9 and 5.10. 
The latter figure is particularly striking, showing the posi- 
tive potentials near the booms distorting the negative poten- 
tials emanating from the aft region. Although the ground 
conductor remains discharged in this case, the solar cells 
on the lower portion of the satellite charged to as much as 
-20 volts, and the reference band charged to -34 volts. These 
results indicate that even in sunlight , SSPM measurements may 
not be accurately referenced to plasma ground. 

Potential contours during discharge of the SCATHA 
model from high negative potentials are shown in Figures 5.11 
through 5.14. Solar illumination was applied using the 
SPINNER model (sun again incident perpendicular to the spin 
axis) after the satellite was charged in eclipse in the high 
temperature double Maxwellian plasma. As the discharge in 
sunlight begins, local regions of strong differential charg- 
ing persist as the overall net charge is dissipated. Some 
surfaces reach positive potentials of up to 2 kV due to the 
transient persistence of differential charging. Figures 5.11 
and 5.12 shew the potential contours after 1 second of sunlit 
discharging, when the persistence of differentially charged 
regions is still clear. After 40 seconds of discharging, as 
■‘llustrated in Figures 5.14 and 5.15, much of the differen- 
tial charging has been removed. The SC6 boom is returning 
from a high positive potential to near plasma ground by col- 
lecting ambient electron current from the plasma, a relatively 
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SCATHA potential contours 
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. Figure 5.10. SCATHA potential contours. 
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SCATHA IN SUNLIGHT; HIGH T SINGLE MAXWELLIAN ENVIRON 
MENT; CONTOUR STEPS = 1 VOLT. 


F:gure 5.11. 


SCATHA potential contours. 
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SCATHA DISCHARGING IN SUNLIGHT; CONTOUR STEPS = 50 V. 


Figure 5.15. SCATHA potential contours. Time = 40 seconds 
after entering sunlight. 
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slow process. This simulation graphically illustrates the 
complicated dynamics involved in the SCATHA exit from eclipse. 
Since many of the interesting processes occur on a timescale 
of less than a full rotation of the satellite, the SPINNER 
model used has undoubtedly distorted the actual time sequences 
somewhat . 


77 


6. SCATHA CHARGING SIMULATIONS USING EXPERIMENTAL DATA 
6.1 DAY 87, 1979 

The SCATHA vehicle charged rapidly shortly after going 
into eclipse on Day 87, 1979, following an injection event. 

The detectors on board SCATHA have transmitted a wealth 
of information on the plasma environment and corresponding 
spacecraft potential during this event. Armed with this in- 
formation, and an accurate representation of the spacecraft, 
we have been able to make the first direct comparison between 
the charging behavior predicted by NASC4V 0 , and that actually 
observed for a real satellite in space. We have also been 
able to show that the physical model upon which NASCAP rests 
is a sound one. 

The environments used in this simulation were those 
described in Chapter 3, and shown in Table 3.1. The repre- 
sentation of the SCATHA satellite used was the "one-grid" 
model described in Chapter 2. 

The potential reached by a spacecraft bathed in a 
plasma environment depends on at least three factors. 

1. The nature of the environment (temperature and 
density) . 

2. The time it has been exposed to the environment. 
(Charging or discharging is not instantaneous.) 

3. The potential of the spacecraft prior to the intro- 
duction of the new environment. 

To properly simulate the response of the spacecraft to the 
charging environment, NASCAP takes all of these factors into 
accour t. After each cycle, the time elapsed is checked, and 
the environment parameters used updated to the most recent 
time for which data was measured. The data points are 
typically 60 seconds apart. 
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The results are shown in Figure 6.1. The NAS CAP 
simulation reproduces the two major jumps in potential , but 
misses the remaining two minor jumps. Quantitative agree* 
ment is excellent considering the sensitivity of the NAS CAP 
predictions to the values of the material properties used. 

The NASCAP simulation is slower to respond to changes in 
environment than the real satellite, because the environment 
changes occur in ^60 second steps rather than the continuous 
adjustment experienced in space. 

In addition, the slow discharge rate predicted, follow- 
ing the two charging pulses, would have been faster if shorter 
computational timesteps had been used. 

The Day 87 simulation is the first real test of both 
NASCAP and the physical model on which it is based. The re- 
markable agreement between the NASCAP predicted potentials and 
those actually observed on a real satellite in an actual space 
environment, shown in Figure 6.1, confirms their validity. 

We can now say with confidence that the physical processes 
which control spacecraft charging are understood. 

6.2 DAY 89, 1979 

The SCATHA SC4-2 electron gun was operated at a variety 
of current-voltage combinations during pass 89-4. These gun 
operations induced a complex response by the spacecraft ground, 
the insulating surface potentials on the SSPM's, and the SC2 
probes. Below we will describe a qualitative picture of the 
beam dynamics, the satellite environment, and the charging 
processes which occurred during pass 89-4. 

6.2.1 Overview and Spacecraft Ground Potential 

The beam dynamics of the electron gun are examined in 

detail in Chapter 7. The electron gun emits monoenergetic 

2 

electrons from an area of approximately 1 cm , at currents 
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Figure 6.1. NASCAP simulated SCATHA charging response for 
Day 87 eclipse. 
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ranging from 0.01 to 13 mA, and voltages ranging from 0.3 to 
3.0 keV. Such a beam is not heavily space charge limited: 
for a 1 kV, 1 mA mode the Child-Langmuir limiting distance is 
8 cm, much larger than the beam diameter. However, the 
escaping beam is significantly spread by its self-field and 
by the satellite field. Electrons returning to the vehicle 
do so after excursions large compared to the satellite radius 
and the beam returns isotropically to a first approximation. 

The relative magnitude of the important charging cur- 
rents is illustrated in Table 6.1. Since the escape of photo 
electrons from the vehicle will be effectively prevented as 
soon as the ground charges positively, the range of currents 
available should be sufficient to charge the vehicle to the 
beam potential in all except perhaps the 0.01 mA case. Ob- 
servations of the actual vehicle potential, as monitored by 
SC10, confirm this prediction. 

TABLE 6.1. SATELLITE ENVIRONMENT 

PHOTOCURRENi’ 

2 

JpH ~ 2 nA/cm 

AMBIENT THERMAL ELECTRON CURRENT 

J_„ - 0.1 nA/cm^ (0_ = 1 keV, n = 1 cm” 3 ) 

In t* 6 

RETURNING BEAM CURRENT 

A SAT ~ 20 m 


*BEAM * 01 


J BEAM ,05 


0.1 


0.5 


1.0 


6.0 


mA 


5.0 30.0 


nA/cm ^ 


The observation that beam currents of 0.1 mA suffice to 
control the satellite ground place limits on the magnitude of 
any low energy electron component in the ambient plasma. As 
the satellite charges positively, the collected low energy 
electron current density for a Maxwell xan distribution is given 
by 



-12 

* 2.7 x 10 



( 6 . 1 ) 


where V is the satellite potential in volts, n the low energy 
density in cm~^ and 0 the low energy temperature in eV. For 

_Q 2 

V * 3 kV, requiring J^ ow << 0.5 * 10 A/cm (see Table 6.1) 
gives 


( 6 . 2 ) 

Thus if n ^ 0.1 cnT^, 0 >> 3 eV, and if 0 % 1 eV, n << 0.06 
cm~^ . Low energy electron components which violate the in- 
equality in Eq. (6.2) would have prevented the beam current 
from charging the vehicle to near the beam voltage. 

6.2.2 Charging Response of Insulating Surfaces 

The surface potential monitors on SCATHA show a com- 
plex response to the range of beam currents employed on 
Pass 89-4. Here we focus on a general description of the 
response of loosely coupled insulating surfaces — these gen- 
eral considerations will be applied to the specific responses 
of the SSPM kapton sample and to the SC2 probe responses in 
Sections 6.2.3 and 6.2.4. 


(-?75 « 6 X 10 



82 


The net electron flux to an insulating surface is 
given by 


J NET " “ J TH “ J BEAM + J SEC + J PH 


( 6 . 3 ) 


where 

* incident ambient flux less backscatter 

JfiEAM * returning beam flux 

J gEC * resulting secondaries (from ambient and 
u beam electrons) 

Jp H * photoelectron flux 

Absolute values are used for the J's with the signs indicating 
the direction of the charging they induce. The term can 

dominate or be comparable to J pH , as indicated in Table 6.1. 
Since J gEC and J pH are characterized by temperatures of a few 
eV, the escape of secondary and photoelectrons is completely 
controlled by the sign of the surface electric fields. Ion 
and conductivity currents are small and can be ignored for the 
purposes of this discussion. 

A qualitative picture of secondary electron production 
must be kept in mind to develop a model of insulating surface 
response. The typical behavior of secondary electron yield 
as a function of incident electron energy is illustrated in 
Figure 6.2. In region I, the surface will tend to charge 
positively. The extent of positive charging is limited by 
the suppression of secondary electron escape as the surface 
field becomes electron attracting. When E RIN > E 2 / the sur- 
face charges negatively until E RIN = E 2 . The value of Ej is 
dependent on the angular distribution of incoming electrons. 

The net flux to a surface will be dominated by dif- 
ferent terms in Eq. (6.3), depending on both the relative 
magnitudes of these terms and on the surface electric fields. 
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Figure 6.3 illustrates the various possibilities for an in- 
sulating surface surrounded by conducting satellite ground. 
Case I occurs when the satellite and the insulator are charged 
to the same positive potential, so that surface fields are 
attracting and secondaries and photoelectrons do not escape. 
The surface then charges negatively with respect to Vg AT . 

The other extreme case. Case IV, occurs when the surface has 
charged towards zero enough to reverse the field on the sur- 
face and cause escape of low energy electrons. Mow the net 
flux is positive and the surface will tend to charge towards 
V SAT* Ne i* her Case I nor Case IV represents a stable equi- 
librium. 


Equilibrium can be reached in two different ways, as 
shown in Cases II and III. If E RIN < E 2 or J BEAM < J pH » then 
equilibrium will be reached as enough low energy electrons 
escape to balance the incident electron currents. Case II 
represents such an equilibrium situation, where V < Vg AT and 
the surface field is reduced enough to limit the required 
fraction of low energy electrons, F. As the surface moves in 
and out of sunlight, J pH will change and the surface poten- 
tial will move correspondingly to adjust F. Thus oscillations 
in surface potential with satellite rotation are expected for 


Case II. If the beam flux dominates J 


PH' 


and if e kin > V 


then limiting of low energy electrons cannot lead to current 
balance. Instead, as in Case III, the surface potential 

E„. Transitions between Case II and 


adjusts until E RIN - ~ 2 . 
Case III are expected wi 
in and out of sunlight. 


Case III are expected when ^ J pH and the surface moves 


Equilibrium in Case II will be reached when the sur- 
face fields are nearly zero, since low energy electrons are 
totally limited with fields of only a few volts per meter. 
Thus the equilibrium voltage results from a purely geometric 
consideration: when does field reversal occur for the in- 

sulator neighboring a conductor? A NASCAP model using a 
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Figure 6.3. Current collection by an insulator backed by 


30 cm rectangular insulating region on a conducting plate 
charged to 3000 volts predicted this reversal to occur at 
500 volts differential. However, it is impossible to obtain 
sufficient resolution using NAS CAP to correctly resolve the 
sharp gradient in charge distribution near the metal/insulator 
interface. A model problem is illustrated in Figure 6.4 
which can be solved approximately to evaluate the field re- 
versal potential. For a conducting annular ring with inner 
and outer radius A and B respectively, field reversal occurs 
when the differential potential at the center, AV, is bounded 
by 


!p(!)< iv <?Ip(!) 

For a 30 cm insulator on a 1 m radius satellite, we find 
150 < AV < 300 volts. 

6.2.3 SSPM Response 

In this section we discuss the response of the large 
kapton sample during pass 89-4. This sample exhibited both 
the largest range and the most consistent pattern of charging 
responses during the gun operations. 

For all but two modes in sunlight, the kapton surface 
stayed within 100 volts of spacecraft ground. These are all 
examples of equilibrium occurring by the Case II mechanism 
described above. 

When the beam energy was 3 kV, then the returning beam 
electrons had greater kinetic energy than the second cross- 
over for their angular distribution. This made the large 
kapton sample charge substantially negative with respect to 
spacecraft ground when the beam current density was greater 
than the photocurrent density. Such was the case for 3 keV 
6 mA beam in sunlight and the 3 keV 0.1 mA beam in eclipse. 
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SSPM voltage profiles for these two examples are shown in 
Figures 6.5 and 6.6. Note the absence of spin period fluctua- 
tions in 6 mA example but their presence# albeit weak# for ^he 
0.1 mA case. The origin of these fluctuations is in magnetic 
field effects on the longer beam excursions in the low current 
case (see Chapter 7) . 

The validity of this surface electric field determining 
current balance model is shown dramatically when the electron 
beam was run in the 3 keV 0.1 mA mode while the satellite was 
in sunlight. In this case when the sample is exposed to the 
sun# the photocurrent is much greater than the returning beam 
current# causing the current balance to be reached as in 
Case II# with the surface nearly at spacecraft ground. How- 
ever# as the sample rotates into shadow# the photocurrent 
disappears and the sample is in a charging beam environment 
and responds accordingly. This strongly spin modulated 
charging is shown in Figure 6.7. 

Not only can one see sunlight effects for the 3 keV 
charging modes, but there is substantial spin modulation in 
the 1.5 keV 0.1 mA and 1.0 mA modes. This is shown in Fig- 
ure 6.8. While the precise mechanism for this modulation is 
not known, two features of the voltage profile are important. 
First# the sample never charges more than 100 volts negative 
with respect to spacecraft ground. This indicates that the 
beam energy is near or below the second crossover, £ 2 # and 
the sample behaves as in Case II. The second, and extremely 
important feature# is that the response is independent of 
beam current, indicating that purely geometrical effects 
determine the electric field structure. This is common to 
the Case II responses; that is, the voltage is not very sen- 
sitive to beam parameters. 
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Figure 6.5. SSPM voltage for 3 kV 6 mA beam in sunlight. 
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6.2.4 SC2 Response 

The response of the SC2-1 and SC2-2 probes to the 
electron beam operations are shown in Figures 6.9 and 6.10. 

The plots stop at the failure during the first 3 keV 6 mA gun 
operation. Certain features stand out dramatically. First 
is the similarity of the response with beam energies of 0.5 
and 1.5 keV for all beam currents. In all these cases the 
spheres were ^350 V negative relative to the spacecraft. 

Indeed there is less than 40 volts total spread for modes 5, 

6, 7 , and 8 observed on SC2-2. This agrees well with a field 
limited behavior as described in Case II , where the potentials 
are determined by geometrical field effects. The 300 volt 
beam, of course, could not produce a greater than 300 volt 
differential. 

When the 3 keV 6 mA mode came on we see that the SC2-2 
charged rapidly 100 volts more negative, then started to drift 
further negative at a rate of about 20 volts per second. This 
type of behavior is typical of differential induced charging: 
that is, where the development of differential charging pro- 
duced saddle point inhibits secondary electrons from escaping 
a surface which would normally be m equilibrium, and as a 
result the surface slowly charges more negatively. This is 
shown in Figure 6.11. (The sudden zero potential reading is 
due to a data dropout.) 

Certain features of the data are not well understood. 

In particular the magnitude of the shadow pulse for the 300 eV 
gun operation is dramatically large compared to the magnitude 
at other times. One possible explanation is that the 100 volt 
incident electrons may be near or below and as such do not 
dominate the dark current collection. Thus the spheres re- 
spond to the weak ambient environment as opposed to the elec- 
tron beam. This is an area which needs to be studied further. 
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6.2.5 Speculation About the SC-2 Failures 


From the SSPM response it is clear that the 3 keV 6 mA 
beam was sufficient to charge kapton to about 1500 volts nega 
tive with respect to satellite ground. Since the booms were 
kapton covered and the spheres were observed at no more than 
500 volts negative, there probably was a 1000 volt differen- 
tial between the booms and the spheres. Consequently, albeit 
a lower differential than associated with surface discharges 
in the laboratory, a kilovolt differential with the appro- 
priate polarity to inject electrons onto the spheres existed 
just prior to their failure. The role of sunlight in trig- 
gering such a discharge is not apparent, but the conclusion 
that differential charging induced discharges led to the SC2 
failures seems warranted. 
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7. ACTIVE CONTROL SIMULATIONS 


The SCATHA spacecraft includes two experiments de- 
signed to achieve active control of the vehicle potential: 
the SC4-1 electron gun and the SC4-2 ion gun. In cons ltation 
with the SC4 experimental team at the Air Force Geophysics 
Laboratory, specific features of the beam operations observed 
during flight testing were identified as of special interest, 
and simulations were performed to elucidate the nature of 
these operations. Sections 7.1 and 7.2 provide estimates of 
the current-voltage characteristics for the SC4-1 electron 
gun, and also assesses the expected magnitude and character- 
istics of the expected return currents. Section 7.3 presents 
an analysis of the observed inability of the electron gun to 
completely discharge the vehicle during a natural charging 
event. Analysis of the operation of the ion gun is complicated 
by the fact that the gun operates in a regime wherein the emis- 
sion is highly space charge limited. The current version of 
NASCAP does not include such effects. A simple model of the 
operation of the ion gun is presented in Section 7.4. 


7.1 SIMULATION OF SCATHA ELECTRON GUN OPERATION 

The electron gun on board SCATHA operates in a regime 
of low current and moderately high voltage. By this we mean 
that the one-dimensional Child's law limiting distance is 
large compared to the bea*v. radius. Since the beam radius 
is only 0.5 cm, the beam dynamics are multi-dimensional as 
demonstrated by the 1-D limiting distance in Table 7.1. Un- 
fortunately, the deviations from one-dimensional behavior make 
any direct analysis intractable. Simulations of the beam 
electrostatic self interactions have shown that the deviation 
from one-dimensional behavior is dramatic, and cannot be 
treated in a perturbative fashion. Indeed, since the self 
forces expand the beam quite rapidly, any similarity to 
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TABLE 7.1. ELECTRON BEAM PARAMETERS AND CORRESPONDING 
1-D CHILD'S LAW LIMITING DISTANCE 


V 

I 

1-D Limiting 

(volts) 

(IRA) 

Distance (cm) 

3000 

6.0 

8 


0.1 

60 


0.01 

190 

1500 

12.0 

3 


one-dimensional behavior is rapidly lost even in the most 
heavily space charge dominated cases. 

For the purposes of analyzing satellite response, we 
must determine certain important features of the beam dynamics. 
The first is the extent of particle excursions. For all cases 
some beam particles join the surrounding plasma. However, the 
variation of satellite potential with beam current indicates 
that, for the higher beam currents, nearly all the beam 
electrons return to the satellite. The data presented in 
Table 7.2 indicates that the ambient plasma current collected 
by the satellite when charged to +3000 volts cannot be larger 
than 0.1 ma. Thus, for the 6 ma beam, over 5.9 ma must re- 
turn to the satellite. The typical orbit length and excur- 
sion time is necessary to estimate magnetic field effects 
and returning beam current density. The mean orbit length is 
determined primarily by the self space charge of the beam for 
all but the lowest currents, where magnetic field effects 
are important. 

Equally as important, but far more difficult to pre- 
dict, is the returning beam density and angular distribution 
as a function of position on the spacecraft. When integrated 
with material surface response functions, this would enable 
a three-dimensional charging analyzer program (e.g., NASCAP) 
to predict the satellite surface potentials. What makes this 

10 J 


TABLE 7.2. SATELLITE POTENTIAL AS A FUNCTION OF CURRENT 
FOR 3 keV ELECTRON BEAM 

1 V 

(ma) sat 

0.01 0 

0.1 3000 

6.0 3000 

difficult is not any subtle nuances of the physics; everything 
can, in principle, be determined using the Lorentz force and 
Poisson's equation. However, even for the 1.5 keV, 13 ma beam, 
the excursions are more than an order of magnitude larger than 
the satellite. Thus, the calculation of the self-consistent 
charge densities and particle orbits with sufficient accuracy 
to predict current variations along the surface is a sub- 
stantial numerical problem. For the longer length orbit, in 
particular for the 3 keV, 6 ma case, the problem is simply 
intractable by straightforward simulation. 

In order to estimate particle orbits, we have devel- 
oped a very simple model which accounts for self space charge 
in a perturbative fashion. The model, which is described in 
detail below, becomes more accurate for larger orbits. To 
demonstrate its validity, we have compared it with a two- 
dimensional simulation of a 1.5 keV, 13 ma beam emitted norm- 
ally from a 1.8 meter radius sphere. (The sphere size was 
chosen to approximate the capacitance of the SCATHA space- 
craft.) 


7.1.1 Simulation of the 1.5 keV, 12 ma Electron Gun 

The 2-D simulation was done using a finite element 
Poisson's code in conjunction with time-dependent particle 
pushing. The satellite potential was fixed at 1.5 keV. The 
mesh outer boundary was at 21.8 meter radius and the boundary 
potential was fixed at zero. While the sensitivity of the 
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calculations to the boundary location was not examined, the 
fact that most trajectories were well within the mesh suggest 
the boundary location was acceptable. The computational mesh 
is shown in Figure 7.1. The calculation was started with no 
space charge in the mesh. As particles were emitted, they 
rapidly developed a space charge field across the beam which 
spread the beam dramatically. This is shown in Figure 7.2. 

The slow electrons at the beam head further enhanced this 
spreading, leading to a major accumulation of space charge 
10 meters out from the satellite. Figure 7.3 shows a super- 
position of all the orbits for t * 0 to 10 usee. Figure 7.4 
shows the electrostatic potentials at that time. 

Although the actual space charge barriers oscillate, 
the calculation was stopped at 15 usee and a single set of 
representative orbits were calculated. These are shown in 
Figure 7.5. They consist of several trajectories correspond- 
ing to different initial positions across the beam width. 
(Since the calculations were performed in R-Z geometry without 
angular momentum about the Z axis, particles could cross the 
Z axis. When this occurs, they appear to be reflected from 
the axis, since a position (-r,z,0°) is equivalent to 
(+r,z,180°) and thus (-r,z) -*■ (r,z) in the R-Z code.) Sev- 
eral important features are seen from those trajectories. 
First, the beam spreads farther than it propagates. This 
implies that the self-expansion forces are dominant in de- 
termining beam dynamics. Secondly, the particles all miss 
the sphere on the first orbit. Two-thirds of the beam finally 
hit the sphere after an orbit of ^450° and ^12 usee. The re- 
maining third appeared to be in a far longer orbit. This 
implies that single pass theories will overestimate the extent 
of excursions. Since only about one percent of the beam can 
actually escape, the potentials and orbits presented here have 
substantial uncertainties and are, at best, qualitative in 
nature. To provide more accurate results would require ex- 
orbitant amounts of computer time. 
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R (m) 

Figure 7.4. Electrostatic potentials after ^10 usee of gun 
operations at 13 ma, 1500 eV. Shown are (A) 
Coulomb potential (upper solid curve) ; (B) Laplace 

potential between 1.8 m and 21.8 m spheres (dashed 
curve) ; (C) potential opposite the beam direction 

(dotted curve) ; and (D) potential in the beam di- 
rection (lower solid curve). The "space charge 
barrier," defined as the maximum difference be- 
tween (B) and (D) , is ^240 volts at a radius of 
\7 meters. 






7.1.2 A Perturbative Model 

As the beam voltage increases and the current decreases , 
the numerical difficulties multiply so rapidly that even the 
simplistic type of calculations presented here become pro- 
hibitive. 

Thus, we are forced to resort to simpler models of the 
beam-space charge interactions. From the current densities in 
the beam, one can see that the strongest forces in the system 
are the 1/r field from the sphere, and the beam expansion 
forces. From Gauss's law these can be seen to be 


(F exp ) 


12x10 


-3 


max 


2iTrVe 


tt*10" 2 2 . 3 X 10 7 8.8*10" 12 


1900 V/m 


(F ) = I 5 ? 0 * 800 volts/meter . 

r max 1 * 8 

The self-retarding fields can be estimated from the numerical 
results (Figure 7.4) by taking the maximum potential difference 
due to the space charge barrier and dividing it by the dis- 
tance to it: 

(Fret 1 : > V V R 1 - v ec (R l )1/R l 


; jZY] 6 ~ 4 5 volts/meter 

By including the first two forces in the particle dynamics and 
approximating the effect of the third by decreasing the initial 
velocity, we have constructed a computational algorithm that 
successfully estimates the extent of particle orbits. Since 
the kinetic energy is not returned to the particles as they 
are deflected from the space charge barrier, the orbits are 
qualitatively incorrect as these particles return to the 
vicinity of the sphere. The radia.. force is estimated via a 
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paraxial theory and is applied only while the particle has 
substantial velocity in the initial beam direction. From the 
orbits of the outermost beam particles, the height of the 
space charge barrier is estimated by integrating the free 
charge Green's function as if the beam were a uniform disk. 
The uniform cross-section approximation is valid only during 
the initial paraxial phase. (Uniform density is a self- 
similar solution to the paraxial beam expansion problem.) 

The calculation is performed iteratively until the height 
of the space charge barrier predicted is equal to the energy 
initially extracted *rom the beam. The resultant electron 
orbit for the 1.5 keV, 12 ma case is shown in Figure 7.6. 
While not in good agreement as the particle returns to the 
sphere, the overall orbit extent is quite good. This pro- 
cedure should improve for higher voltage, lower current 
beams where 


F, . >> F >> F , 

1/r exp ret 


The predicted transit time was 8 microseconds. The transit 
time is important in determining the effect of magnetic fields. 
For pass 89-4, 1979, the approximate field strength was 100y 
or 10" 3 Gauss. This corresponds to an electron w c of 

pQ 4 

uj = — - 2 x 10 rad/sec . 
ce me 


Consequently, the effect of the magnetic field on the orbits 
is negligible, co ce T = 0.2 radians at most. 

We have applied this model to cases of 3 keV, 6 ma and 
3 keV, 0.1 ma electron beams. Resultant orbits are plotted 
in Figures 7.7 and 7.8. The 6 ma beam was limited by its own 
space charge some 50 meters from the satellite. The particle 
transit time is approximately 30 microseconds. For these 
particles u> t is 0.6 radians which is only enough to make 
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z-nxis 

Figure 7.6. Approximate electron orbit for 1500 eV, 12 ma case 




Approximate electron trajectories for 3 keV, 0.1 raa 
case in the absence of a magnetic field. 


very small deviations in the orbit. Figure 7.9 shows the 
trajectories assuming B is normal to the plane of the orbit. 
The effects are barely discernible. However, for the 0.1 ma 
beam, the magnetic field dominates. Indeed, in the low cur- 
rent limit the beam is magnetically limited for most orien- 
tations. Figure 7.10 shows an example of a magnetically 
limited orbit. 

7.1.3 Conclusions 

The computer results suggest that, to a reasonable ap- 
proximation for all but the lowest current cases, the beam 
electrons return to the sphere uniformly. If, in addition, 
we make the further simplification that the outgoing beam is 

spherically uniform, we can use the Langmuir-Blodgett theory 

[231 

of a spherical diode 1 J to make a simple estimate of the ex- 
cursion distance and time of the beam electrons. 

Under the above assumptions, the beam current (approx- 
imated as 100 percent returning) , beam energy, and travel 
distance are related by 

2i = 2.93 x 10” 5 V 3/2 j (-a) 2 

where the factor of 2 accounts for the beam's contributing 

2 

space charge both leaving and returning, and (-a) is a 

function of the ratio of spacecraft radius, r , to beam ex- 

s 

cursion radius, r B , given tabular ly by Langmuir and Blodgett. 
The mean bean excursion time. At, is found by requiring the 
total beam space charge to equal the spacecraft charge: 

i4t * 4,ie o v (s- - ib) • 

' S B ' 
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Les for 3 keV 





Figure 7.10. Approximate electron trajectories for 3 keV, 0.1 ma case in 
the presence of a magnetic field. 



Table 7 . 3 presents the excursion distance and time for 
several cases calculated by all three methods. We see that 
the perturbative method, described in the previous section, 
compares quite well with the simulation and should be an ex- 
cellent estimate for the low current high voltage beams. How- 
ever, the Langmuir-Blodgett based estimates are good to a 
factor of two and are a ready method for predicting excursion 
times. 


TABLE 7.3. EXCURSION DISTANCES FROM A CONDUCTING SPHERE 
CALCULATED VIA SIMULATION, PERTURBATIVELY AND 
USING THE LANGMUIR-BLODGETT ESTIMATE 


V 

I (ma) 

R max lR - Z) 
Simulation 

R Pertur- 

max 

bation 

max 

(L-B) 

At (usee 

1500 

12 

17 

18 

29 

25 

3000 

6 


44 

90 

100 


1 



290 

600 


.1 


540 

■vlOOO 

0/6000 

300 

.1 



135 

150 


Comparing these excursion times to the gyroperiod 2tt/cj 'v 

G 

300 usee, we see that space charge effects dominate magnetic 
field effects for currents above 1 ma for a 3 kV beam, and 
currents as low as 0.1 ma for a 300 volt beam. Interestingly, 
only for the 3 kV, 0.1 ma case were magnetic field effects 
visible in the SSPM response. 
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7.2 


CURRENT-VOLTAGE CHARACTERISTICS OF THE SCATHA SPACE- 
CRAFT 


Ambient current levels versus satellite potential have 
been obtained for a simple SCATHA model. The results are very 
sensitive to the environmental description employed. Two 
cases were considered which should bracket the actual environ- 
ment: (1) a low energy dense plasma, described by single ion 
and electron Maxwellian distributions with n fi * n^ • 10 cm 
T e “ T i * 2 eV; an< * (2) a high energy plasma, also Maxwellian, 

with n ■ n. ■ 0.2 cm”^, r ■ 8 kV, T, » 13 kV. NASCAP was 
ex e x 

used to calculate ambient current to the four-grid model of 
the satellite assuming uniform charging in eclipse. The cur- 
rent versus voltage characteristics are shown in Figure 7.11. 
The area between the two curves represents the range of pos- 
sible I-V characteristics expected. As discussed above, when 
the beam current reaches the level of the ambient current, the 
satellite potential will be controlled at beam potential and 
the return current will increase sharply. The curve for the 
low energy plasma certainly overestimates the ambient current 
due to the probe model employed for collection. It is there- 
fore encouraging that the beam current available, 1 yA to 
13 mA, includes the spectrum of expected I-V characteristics. 

The results in figure 7.11 were obtained for the satel- 
lite in eclipse. In sunlight, the ambient current is domi- 
nated by photoemission. NASCAP was yjed to predict the net 
current in sunlight, and the results vary between 60 yA for 
a uniformly negatively charged (or discharged) satellite to 
18 uA for a differentially charged satellite (with the con- 
ductors held negative) . Photoemission is of course completely 
limited when the satellite charges to positive potential. 
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7.3 POTENTIAL BARRIER FORMATION ABOVE THE SC4-1 GUN 

When SCATHA encountered a magnetospheric substorm 
which charged the vehicle to -6 kV, an attempt was made to 
reduce the charging using the SC4-1 eiactron gun operating 
at 150 volts and 1 mA. The attempt was only partially suc- 
cessful, in that it was possible to drivs the vehicle poten- 
tial only to -1 kV, but not higher. One reason for the in- 
ability of the gun to completely discharge the satellite 
could be the formation of a potential barrier above the 
electron gun which would prevent the emitted electrons from 
escaping. In a strongly charging environment, the ground 
conductor at -1 kV around the gun aperture will be a region 
of rexatively high potential, so the formation of a barrier 
is at least plausible. We have used the NASCAP code in an 
analysis of three questions regarding this hypothesis: 

1. Will such a barrier form? 

2. How much time is required to form a barrier? 

3. Where will the barrier cause returned particles 
to hit the vehicle? 

We began a simulation using the four-grid SCATHA model 
which had been charged in eclipse using a high temperature 
double Maxwellian environment described in Chapter 4; the 
ground potential was -6410 volts. We then fixed the ground 
conductor to -1000 volts while maintaining any existing dif- 
ferential charging* the resulting potential contours are 
shown in Figure 7.12. The SC4-1 electron gun is located 
on a conducting region of the bellyband, so that ejected 
electrons are emitted at -1 kV less the gun potential. 
Initially, the external fields above this region are repel- 
ling, so that all 150 volt electrons could escape. The ex- 
posed dielectrics were well above their equilibrium potentials 
in this environment since they followed the +5 kV jump in the 
ground conductor. The ground conductor was then held at 
-1 kV while the exposed dielectrics and other conductors 
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differentially charged. A barrier above the SC4-1 gun then 
slowly formed, chiefly as a result of the charging of both 
the neighboring solar cell covers and the booms. 

The potential contours after 395 seconds of charging 
are shown in Figure 7.13, and the associated trajectories of 
emitted electrons are shown in Figures 7.14 and 7.15. A barrier 
of approximately 141 volts existed at this stage, enough to 
deflect the trajectories significantly but not to cause any 
particles to return. After 415 seconds, the barrier height 
increased to 144 volts. Now 40 percent of the emitted 
electrons returned to the vehicle, as shown in Figures 7.16 
and 7.17. Notice that many particles escape by glancing off 
the barrier, since the barrier height diminishes rapidly away 
from the bellyband. By 435 seconds, the barrier increased 
to about 150 volts, and 80 percent of the particles returned 
to the vehicle, as shown in Figures 7.18 and 7.19. 

Emitted electrons which return travel along the belly- 
band, away from the region of high negative potential gen- 
erated by the SCI teflon mask. These electrons cannot pass 
the potential barrier around the SC6 boom, and they return to 
the vehicle on the bellyband between the SC4-1 gun and the 
SC6 boom or on the solar cells above this region. 

The total time for the barrier to form in the above 
simulation was 8 minutes. However, this time is very sensi- 
tive to the initial charging conditions, which is in turn 
sensitive to the material properties, and these are poorly 
known. The "switching time" for a change from under 25 per- 
cent current returned to over 75 percent returned is about 
30 seconds. Both the above times scale inversely with the 
assumed ambient density, which was 0.2 cm 3 in the above case. 
Space charge limiting of the emitted beam (neglected above) 
could also allow more rapid barrier formation due to local 
charging of solar cells near the SC4-1 gun. 
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Figure 7.13. Potential contours after 395 seconds. Note the 
potential barrier above the SC4-1 gun location. 
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Figure 7.14. Trajectories of 150 volt electrons from SC4-1 
gun after 395 seconds. All particles escape. 
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Figure 7.18. 


Trajectories of 150 volt electrons from SC4-1 
gun after 435 seconds. Eighty percent of emitted 
particles return. 
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7.4 


SPACE CHARGE LIMITED ION BEAM EMISSION 


One of the experiments on board the SCATHA satellite 
is a plasma discharge ion emitter. The working gas is Xenon , 
whose large mass makes the unneutralized beam extremely space 
charge limited at the operating voltages and currents. In 
this section we present an analysis of the emitter based upon 
simple space charge limited diode theory. While the model 
developed has a large number of limitations , the satellite 
voltages it predicts are in rough agreement with experiment. 
The model identifies the relevant physical mechanisms of beam 
emission and thus is useful in planning more elaborate multi- 
dimensional calculations. 

7.4.1 Model Description 

We assume the emitter to be a 1.27 cm diameter disk 
held at spacecraft ground. The beam is assumed to be mono- 
energetic with the ion velocities all normal to the emitting 
surface. Ion energy is 1100 eV. At the three current levels 
evaluated (2 mA, 0.3 mA, 0.08 mA) the limiting distances are 
much shorter than the assumed plasma Debye length of 700 cm. 
This length corresponds to a 0 = 1 eV, n g * 1 era 3 background 
plasma. As will be shown, the results are not sensitive to 
the estimates of the background plasma. 

The severe limiting is used to divide space into two 
regions. The first region is between the emitter and the 
limiting layer. The second is outward from the limiting 
layer, which is a virtual anode for ion emission into the 
plasma. In Figure 7.20 we illustrate qualitatively the be- 
havior of the potential along a line normal to the emitter 
surface. 

In Region 1 the ion beam forms a potential barrier 
which reflects the beam back to the emitting surface. Since 
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Figure 7.20. Qualitative dependence of the potential along a 

line normal to the emitting surface. $$at t ^ e 
satellite ground potential, is the virtual 
anode potential, V b is the energy of emitted ions 
and X T is the distance to the virtual anode. 





for all the cases studied here, the satellite current is much 
less than the beam current, the barrier must be sufficient to 
limit the beam. Using our assumption of a monoenergetic beam 
we get 


^A " ^SAT 


where is the beam voltage, $ A is the virtual anode poten- 
tial and 4 > sat is the satellite ground potential. Since the 
initial beam ion density is many orders of magnitude larger 
than the ambient electron density (neutralizer off) , in the 
analysis below we will neglect the effects of ambient screen- 
ing in Region 1. 


In Region 2 we have cold ions leaving the barrier and 
streaming into space. Note that the ion velocities mono- 
tonically increase with distance beyond the virtual anode, 
located at a distance X L from the emitter. The potential 
decays at long distances due to the ambient plasma screening. 
While the screening is essential for any ion current to 
exist, we will show below that the amount of current is in- 
sensitive to both the plasma Debye length and the precise 
nature of the shielding. 


The two regions are coupled by solving for the satel- 
lite potential, $ SAT * for which the ion current in Region 2, 


balances the net satellite current I 


SAT* 




+ I SAT (i SAT ) 


0 


I is gotten from the I-V characteristic of the satellite 

iD X 

as calculated by NASCAP. 
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REGION 1 


The formation of the virtual anode in front of the 
disk emitter is inherently a two-dimensional problem, since 
the space charge expands the beam as well as slows it. For 
the purposes of this study we have developed a simple, albeit 
ad hoc, model of the beam spreading which allows us to esti- 
mate the size and shape of the virtual anode using one- 
dimensional spherical diode theory. We emphasize that this 
is probably the weakest part of the theory, but qualitatively 
it does account for the beam spreading. 

In Figure 7.21a we illustrate the virtual anode forma- 
tion. In Figure 7.21b we show the geometric construction 
which we use to convert the actual multidimensional case into 
a solid angle of 2it(1-cos9) of a concentric spherical emitter. 
The chord of the virtual anode is set equal to the emitting 
disk diameter plus twice the limiting distance. From this 
assumption we get the following geometrical relations: 


r 0 * r l + X L 


r^ = a + 


2X. 


i - cose = — 

* A 


where r Q is the virtual anode radius, r^ the emitter radius, 
X. the limiting distance, a the disk radius, and e the half 
cone angle. We complete the equations by relating the beam 
current I B to the total Child's law current, I g , of a concen- 
tric spherical diode with ratio of radii rg/r^ 


X S = 


2.9 


10 


-5 V 


3/2 


VI > 2 (^) 
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Figure 7.21a. Virtual anode formation in front of an emit- 
ting disk. 



Emitting Surface 


Figure 7.21b. 


Our model emitter where the disk now represents 
a solid angle 2tt(1-cos9) of a sphere of radius 
r^ emitting and forming a virtual anode at a 
radius r-> = ri + X L . 
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T , 2tt (1-cosE) _ 

REG 1 4rf A S 

I PEG 2 " 2 X B I SAT 

The current in Region 1, is almost twice the beam cur- 

rent, I 0 , because most of the current is reflected from the 
space charge barrier. This is especially true in the 2 mA 
case since the beam current, 2 mA, is much greater than the 
satellite current ^0.02 mA. 

The solution of these equations for ■ 1100 volts, 
a * 0.635 cm and M ■ 13.1 amu is given in Table 7.4. Notice 
that only for the high current case is the limiting distance 
smaller than the disk radius, a. The virtual anode becomes 
more hemispherical with decreasing beam current. For the 
0.08 mA case the half cone angle is almost 80 degrees. 

TABLE 7.4. VIRTUAL ANODE PARAMETERS AS A FUNCTION OF BEAM 

CURRENT FOR V. = 1100 VOLTS: EMITTER DISK RADIUS 
a = 0.635 cm 


I B (ma) 

X L (cm) 

2.0 

0.30 

0.3 

1.08 

0.08 

2.92 


R 0 (cm) 

( 1-cosG ) 

1.57 

0.19 

1.89 

0.57 

3.62 

0.81 
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REGION 2 


In Region 2 we have calculated the current emitted by 
a sphere into a plasma. In the cases of interest X Q >> r Q , 
the effective radius of the emitter. Poisson's equation in- 
cluding both beam and ambient plasma densities is 


- I 1 » ghprical . 4ir 

rJ - ♦) 

We approximate the ambient charge 
term 


p ambient 


density using the linear 


** ^ambient 



We have also used the approximation 


-4it p 


ambient 


X D 1 + 




371 


to account for electron acceleration in the sheath. 

The current for a given t > ^ that satisfies the boundary 
conditions at infinity -► 0) is remarkably insensitive to 
the form of ambient plasma screening. Indeed for X D /r Q ■ 200 
the two forms predict the same current to within 1 percent. 
This is due primarily to the dominance of beam space charge 
at small radii and the insensitivity of concentric spherical 
emitters to the radius of the virtual anode. Table 7.5 shows 
the space charge limited current for both approximations as 
a function of X n /r for a spherical emitter at 700 volts. 

uo 3 /2 

The results scale as V and depend only on \ D /r Q , the same 
as a conventional spherical diode. For planar geometries, 
the relationship between the screening length \ D , the emitter 
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voltage and the current density can be solved analytically 
(see Appendix B) . Then the analogy between X D and the asso- 
ciated diode spacing is exact with the effective gap distance, 

Effective' just a con8tant times X D 

2/2 

^effective * 3 ^0 

2 

V effective ^A ^o * 3 ^A 


J 


planar 


2.64 x 10 



with V in volts and d in cm. 

The spherical current is related to the satellite cur- 
rent by the subtended solid angle: 

I SPHERICAL * I SAT /(1 ‘ COS0) 

Examining Table 7.5 it is apparent that the solid angle of the 
emitter is more important than the emitter radius in deter- 
mining the voltage required for a given current emission. 

This is because the largest space charge effects occur close 
to the emitter surface. 

TABLE 7.5. TOTAL SPHERICAL XENON ION CURRENTS INTO A PLASMA 
FOR AN EMITTER AT 700 VOLTS. 


\ /r 

D x 0 

v D ' (mA) 

'(,2 177372 ) 

' D v ' (mA) 

1 

2.59 

1.78 

5 

0.72 

0.64 

10 

0.52 

0.49 

200 

0.25 

0.25 

700 

0.21 

0.21 

Note that at large '• D /r c the current is msensiti 
to both the ratio of the Debye length to emitter 
radius and to the form of the plasma screening. 
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7.4.2 Results and Discussion 

We couple Regions 1 and 2 in order to determine the 
virtual anode voltage required for emitters with radii and 
solid angles as described in Table 7.4 to emit prescribed ion 
currents into a 700 cm Debye length plasma. The satellite 
voltage is found assuming it is 1100 volts below the virtual 
anode potential. The results for various satellite currents 
are shown in Table 7.6. The blank entry for case 1, I gAT * 

0.05 mA is because analysis suggests this is a physically 
unrealizable case. The predicted satellite voltage is posi- 
tive which would severely reduce the satellite current below 
0.05 mA. 

# 

The range of ambient currents to the SCATHA satellite 
in sunlight was determined using the NASCAP code. A mild en- 
vironment was considered, described by Maxwellian electron 

and ion distributions with n^ = n. = 1 cm ^ and T = T. = 1 eV. 

e i e i 

In such an environment, the ambient current is dominated by 
photoemission. The maximum current is emitted when the satel- 
lite is uniformly negatively charged. The current is +60 yA 
at -1000 volts, decreasing slowly with voltage to approxi- 
mately 50 yA at -100 volts. Lower current levels are obtained 
when the insulating surfaces are discharged so that positive 
fields develop to limit photoemission. This occurs in several 
seconds. With the conductors fixed to -1000 volts and the 
insulating surfaces largely discharged, the net current was 
reduced to 18 yA. 
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TABLE 7.6. SATELLITE VOLTAGES FOR THREE BEAM CURRENT LEVELS 
AS A FUNCTION OF SATELLITE CURRENT 


Case 

JL 

' B ■ 2 ■*' V B 

'sat (mA) 

= 1100 Volts 
4 > a (volts) 

^AT (volts) 



0.05 

1288 

- 



0.02 

696 

-402 



0.01 

440 

-660 

Case 

_2 

I_ = 0.3 mA, 

O 

V B = 1100 Volts 




'sat <"*> 

4» a (volts) 

*SAT lvoltsl 



0.05 

570 

-530 



0.02 

30 3 

-791 



0.01 

125 

-905 

Case 

_3 

** 

00 

o 

• 

o 

II 

CO 

HH 

V_ = 1100 Volts 

o 



1 SAT (mA) 

4> a (volts) 

"'SAT <volts) 

0.05 

440 

-660 

0.02 

240 

-860 

0.01 

150 

-950 
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7.5 


DISCUSSION 


The analysis presented here is qualitative in nature 
and was performed for the purpose of identifying the relevant 
physics of space charge limited ion emission from a satellite. 
The results agree with experiment to the extent that lower 
beam currents produce larger negative satellite potentials. 
Also, the magnitude of the potential variations are similar 
to those observed. Our analysis indicates that the magnitude 
of the potential variations is a function of the geometry of 
space charge barrier, and since we estimate this rather 
crudely, our results are not to be considered qualitatively 
accurate. 

More accurate results could be obtained in the follow- 
ing manner. First a local two-dimensional particle pushing 
analysis of the space charge limiting of the emitter would 
be performed. Since the plasma boundary conditions influence 
this only very weakly, this part could be easily performed 
including such phenomena as gun optics. Then given the 
space charge barrier, multidimensional space charge limited 
currents in space beyond the barrier could be calculated as 
a function of satellite potential. This would include the 
dipolar nature of the space potential in Region 2, which will 
increase the beam divergence substantially in the high beam 
current cases. 
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8 . DETECTOR MODELING 


The SCATHA spacecraft has on board several particle 
detector experiments which monitor particle fluxes at geo- 
synchronous orbit. The experiments SC2, SC5 , SC6 , SC7, and 
SC9 were designed to measure electron and ion fluxes in the 
energy range which causes the charging of exposed spacecraft 
materials, 0 to ^80 keV. The SC6 experiment never operated 
in orbit, and the SC7 experiment operated only briefly, so 
attention will focus on the remaining three experiments. 

This chapter describes the use of the NASCAP "DETECTOR" mode 
to simulate the response of the SCATHA particle detectors. 

A particle detector with an unobstructed field of view 
operating on an uncharged, magnetically clean spacecraft can 
measure particle fluxes which accurately reflect the state 
of the ambient plasma environment. Differential charging 
and local geometric effects produce electrostatic fields 
near a spacecraft which can distort the trajectories of in- 
coming particles. Low energy electrons emitted from highly 
charged regions of a spacecraft can be observed as high 
energy electrons incident elsewhere on the vehicle. The 
NASCAP DETECTOR model was designed to assist in the identifi- 
cation and interpretation of such occurrences. The operation 

f 2 1 

of the DETECTOR routines is described in detail elsewhere. 
Briefly, numerical reverse trajectory tracking is used to 
connect particle orbits incident at a specified detector 
location on the vehicle with a postulated phase space distri- 
bution at a large distance. The predicted response can be 
displayed as a function of view direction or incident parti- 
cle energy, and plots of the trajectories are produced. 

The simulations described in the following sections 
were performed to illustrate expected detector responses 
using the charged states of the SCATHA spacecraft generated 
in the studies described in Chapter 5. Three representative 
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cases were selected: 


Case Is High temperature plasma in sunlight. See 
Section 5.3, Figures 5.9-5.11. Spacecraft 
ground at +4 volts. 

Case 2: High temperature plasr.-* in eclipse. See 

Section 5.1, Figures 5.3 and 5.4. Spacecraft 
ground at -6400 volts. 

Case 3: Moderate temperature plasma in eclipse. See 

Section 5.2, Figures 5. 6-5. 8. Spacecraft 
ground at -140 volts. 

The first case represents a nearly uncharged state of the 
vehicle, while the next two cases had significant differen- 
tial charging, especially along the SC2 and SC6 booms. De- 
pending on the operational state of the satellite (i.e., 
whether or not the boom mounted experiments were grounded, 
biased, or allowed to float at the time of inteiest) the 
differential charging pattern represented by these cases can 
be altered significantly. The simulations described in the 
following sections were chosen to illustrate the most severe 
type of vehicle perturbations which can be expected. In all 
three cases, the phase space distribution at infinity was 
assumed to be isotropic and characterized by the appropriate 
single Maxwellian representation from Table 4.1. A constant 
B-field of 10 ^ gauss along the +y axis was used in all cases 
as representative of geosynchronous conditions. Since even 
a 10 volt electron has a Larmor radius of 12 m for this field, 
the turning of orbits by the magnetic field is negligible. 

For further discussion of this point, see Section 8.2 below. 

8.1 SC2 DETECTOR SIMULATIONS 

This section describes simulation of the response of 
the SC2-3 detector mounted on the SCATHA bellvband, near the 
base of the SC2-2 boom. Two similar detectors are mounted 
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at the ends of the SC2 booms: the SC2-2 detector has an un- 

obstructed field of view, while the SC2-1 detector looks back 
towards the vehicle body. As the results below demonstrate, 
the electrostatic fields from the SC2-2 boom can significantly 
perturb the trajectories of particles observed by SC2-3. 

The predicted response of SC2-3 in the Case 1 environ- 
ment is illustrated in Figure 8.1, and the associated incoming 
electron and proton trajectories are shown in Figure 8.2. Since 
the satellite is almost uncharged, the trajectories are nearly 
unaffected. Only the lowest energy electrons (a few eV) are 
significantly perturbed as they travel nearly parallel to 
the SC2 boom. Figure 8.3 displays expected response for the 
Case 3 environment. Since the vehicle ground is charged to 
-6400 V, protons below this energy are not observed. Any 
incident protons below this energy must have originated 
elsewhere on the vehicle, as shown in Figure 8.4. Electron 
trajectories are distorted only slightly in the repelling 
fields for this case. Displayed in Figures 8.5 and 8.6 are 
corresponding results for the Case 3 environment. The SC2 
boom tips are highly charged in this case, while the ground 
conductor is at -l^O volts. Both electron and ion trajectories 
are severely distorted in the spacecraft fields, and the ex- 
pected detector response is accordingly modified. Some low 
energy protons actually orbit the boom tip before arriving 
at the detector, completely distorting the information re- 
garding angular dependence of incoming particle trajectories. 
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6. Trajectories of electrons (top) and protons (bottom) 
logarithmically spaced from 10-19,000 eV, observed 
at SC2-3 location for Case 3 environment. Many pro- 
ton trajectories orbit the negatively charged SC2 
boom tip without actually striking the boom. 
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SC 5 DETECTOR SIMULATIONS 


The SC5 experiment includes two electrostatic analyzers 
which measure electron and ion fluxes between 50 eV and 60 keV. 
These detectors view parallel and perpendicular to the satel- 
lite spin axis, with the bellyband (perpendicular) detector 
rotating approximately in the plane of the ambient magnetic 
field. During magnetic substorms, the response of the SC 5 
bellyband detectors on SCATHA have exhibited variation depend- 
ing strongly on the azimuth about the spin axis. The con- 
siderations presented bel yw agree that such a feature of the 
response is characteristic o' anisotropy of the particle flux 
beyond the range of influenr. R of the satellite and is not 
induced by complex interaction involving the electric fields 
of the charged satellite and the ambient magnetic fields. 

More precisely, the satellite electric fields will not cause 
an isotropic distribution of particles originating from be- 
yond R to appear as an anisotropic distribution at the 
detector. 

Consider first the relevant length scales of the 

problem 


R - X D = 


cm Debye length 

(8.1) 

A « 1.7 
P 

f e 

- i° 2 Vr" 

cm proton Larmor radius 

(8.2) 

3.7 

* 

B Cm 

electron Larmor radius 

(8.3) 


where d eV is the plasma temperature, n g the electron den- 
sity, B the ambient magnetic field in gauss, and Ep(E e ) the 
energy in eV of protons (electrons) in the plasma beyond R. 

For 6 -v 10 4 eV, B ~ 10' 3 gauss, E ^ 100 eV, r. 'u 1 cm -3 

P 
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(6.4) 


X D << 


X D 


* 4 x IQ 


-3 


Even for electrons with sufficient energy to reach the satel- 
lite (E^ 0 % ! v sat |) 


X D /X e * 0,2 (6.5) 

Consequently, the trajectories of protons after entering the 
range of influence of the satellite are negligibly influenced 
by the ambient magnetic field. We therefore ignore B in sub- 
sequent considerations. 

Now consider the eclipsed satellite in a substorm 
environment. Examination of various time scales for dif- 
ferential charging as well as detailed NASCAP calculations 
indicate that the electric field structure in the space 
around the spacecraft when viewed in the frame of reference 
rotating with the stellite, does not change substantially 
in the course of one rotation. A gain in flux that is iso- 
tropic beyond R would appear isotropic at the detector, and 
any apparent anisotropy must be present in the particle 
distributions at large distances from the satellite. 

A primary effect of the satellite induced electric 
field is to distort the volume elements in phase space. To 
determine the region of velocity space beyond R that cor- 
responds to regions sampled at the detector position can be 
determined by reverse tracking of particles entering the 
detector. The trajectories presented in Section 8.4 for the 
SC9 high energy detector are very similar to the corresponding 
simulation results for the parallel SC5 detector, except that 
the very low energy range (-'50 eV) is not monitored by SC5. 
Figures 8.7 through 0.10 illustrate simulated response for 
the SC5 bellyband detector in the Case 2 and Case 3 environ- 
ments. (The trajectories for the Case 1 environment are 
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essentially undistorted straight line paths.) Figure 8.7 
ohows the characteristic cutoff of protons below the vehicle 
potential for Case 2 , and the perturbation of the proton 
trajectories by the charged satellite is visible in Figure 
8.8. In Case 3/ the charged SC2 boom tip ha3 similar effect 
on the incident protons and electrons as was observed for 
the SC2-3 detector. The simulated response and the cor- 
responding trajectories are shown in Figures 8.9 and 8.10. 
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8.3 


SC 7 DETECTOR SIMULATIONS 


The SC7 light ion mass spectrometer measures the cold 
(0 < 100 eV) ionic component of the plasma using three sen- 
sors located parallel, antiparallel, and perpendicular to the 
satellite spin axis. The SC7 detectors have a much larger 
field of view, 110°, than the other particle spectrometers, 
o-i-nce low energy ions are observed, their trajectories can 
be perturbed by the electrostatic fields of the satellite 
even when the vehicle is not highly charged. Two comparison 
cases have been run to illustrate this effect. These tra- 
jectories differ from others presented in this chapter in 
that the energy is fixed for all trajectories while the azi- 
muthal angle of the incident particles is varied from 0 to 
360°. The polar incident angle was fixed at 45°. Fiqure 
8.11 shows trajectories of 50 eV protons entering the SC7-1 
bellyband detector for the Case 1 potentials. Little per- 
turbation of the trajectories is observed in these weak 
fields. Trajectories for the Case 3 moderate temperature 
plasma are illustrated in Figure 8.12. Here significant dis- 
tortion of incoming trajectories results from the influence 
of the highly charged SC2 boom tips. Similar effects can be 
seen for the SC7-2 and -3 detectors. 

8.4 SC9 DETECTOR SIMULATIONS 

The SC9 UCSD charged particle experiment consists of 
three detectors mounted together on the forward surface of 
the satellite. There are two rotating detectors for electrons 
and ions, one scanning the range 0.2 eV to 1550 eV and the 
other scanning the range 1 eV to 81 keV. The third detector 
scans the low energy range for ions only. This section 
reports simulations of detector response at fixed detector 
angles while scanning the lower and higher energy ranges. 
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Detector response in the two high temperature plasma 
cases are represented in Figures 8.13 through 8.16. For 
Case 1, with the vehicle essentially uncharged, the incident 
spectrum to the high energy detector looking parallel tc the 
spin axis appears as in Figure 8.13. Only the very lowest 
energy electron trajectories, at a few eV, are perturbed by 
the weak spacecraft fields, as shown in Figure 8.14. Cor- 
responding results for the highly charged eclipse state. 

Case 2, are shown in Figures 8.15 and 8.16. Protons below 
the spacecraft potential, -6400 V, have trajectories termi- 
nating elsewhere on the vehicle. Finally, the response to 
the moderate temperature environment, Case 3, is illustrated 
in Figures 8.17 and 8.18. The proton spectrum cuts off at 
the spacecraft potential, -140 volts, and the trajectories 
are dramatically perturbed by the SC2 boom. The smooth 
spectrum for incident electrons gives no hint that the low 
energy trajectories are influenced by stray fields from the 
OMNI antenna. 

The response of the low energy ion fixed angle de- 
tector is even more interesting in the Case 3 environment. 
The predicted incident spectrum has a gap in the middle 
energy range, shown in Figure P.19. The corresponding tra- 
jectories, Figure 8.20, illustrate the reason for the com- 
plex response. High energy protons travel to the detector 
in nearly straight paths. Middle range protons are accel- 
erated towards the charged SC2 boom, where they impact, 
leading to the gap in the spectrum. The lowest energy 
electrons are accelerated so strongly towards the boom tip 
that they pass below the boom plane and then orbit the boom 
once without striking, then pass away from the boom region 
and connect tc trajectories at infinity. 
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Figure 8.15. 

Simulated response of SC9 NS detector to incoming 
electrons and protons, 1-81,000 eV, in Case 2 


environment. 
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Figure 8.17. Simulated response of SC9 NS detector to incoming 
electrons and protons, 1-81,000 eV, in Case 3 
environment . 
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9. HOT SPOT THEORY 


The study of dark side potentials in limiting the ef- 
fectiveness of photoelectrons to discharge sunlit surfaces 
has raised the question as to the effect of small shaded in- 
sulators on nearby photoemission. Some estimates of such 
effects have been found by numerical calculations, but in 
order to get a better overall understanding we have also ex- 
amined the problem analytically. The analysis described below 
should prove useful in understanding the importance of highly 
charged hot spots on the SCATHA spacecraft. 

We have constructed a model problem of a single circu- 
lar region on an infinite insulating plane. The circle is 
uniformly charged to a large negative voltage V, (^-5 kV) 
while the remainder of the plane is exposed to intense sun- 
light and an incident high energy electron flux (e.g., 5 kV 
Maxwellian) . In the absence of the negative region, the photo- 
emission would completely discharge the insulator and the 
potential would be a few volts positive uniformly on the plane. 
The "hot spot" sets up a local potential barrier which prevents 
photo and secondary electrons from leaving the insulator. As 
a result, charge accumulates and the insulator goes negative. 
Below we discuss using analytical models of the spatial de- 
pendence of the potential caused by the "hot spot" effect. 

9.1 HOT SPOT IN THE ABSENCE OF SHEATH CONDUCTION 

In the case where sheath conduction is small, the in- 
fluence of a hot spot is similar to that of the negatively 

charged dark side in the sunlit sphere case discussed by 
f 2 1 

Katz, et^ al . The essential feature is that pbotoemission 

and secondary emission can support only very small electric 
fields normal to the emitting surfaces. As a result, equi- 
potential contours from a flat circular hot spot on an infi- 
nite insulating plane must be approximately normal to the 
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plane. The net effect is to make the insulator "disappear" 
(see Figure 9.1). A similar effect with regards to photo- 
emitting booms was observed by Schnuelle, et al. on the NASCAP 
model of SCATHA. 



Figure 9.1. Potential contours around an isolated disk. 

Since the contours are perpendicular to the 
dashed line, the electric field normal to the 
line is precisely zero. 

The analytical potential in the plane of a disk at 
potential V Q of radius a falls off in r as 

2V 

V(r) = arcsin (a/r) (r > a) (9.1) 

7T 

In the limit of r >> r^ (see Table 9.1) 

(9.2) 

From previous experience with photosheath behavior we 
imagine this to be a good representation of the hot spot i' - 
fluence, and a rigorous upper bound of the potential devia- 
tion from the same case without a hot spot. 



TABLE 9.1. COMPARISON OF ANALYTICAL REPRESENTATIONS OF 
HOT SPOT POTENTIALS 


r 

2/ir arcsin 1/r 

2/irr 

1 

1 


1.001 

0.991 

0.637 

1.01 

0.910 


1.1 

0.726 


1.2 

0.627 

0.531 

1.5 

0.465 

0.424 

2.0 

0.333 

0.318 

3.0 

0.216 

0.212 


As a check of the usefulness of this approach we set 
up a NASCAP model of a charged square on an insulating octa- 
gonal surface. Increased secondary electron emission was 
used instead of sunlight so that both sides of the octagonal 
object would be the same. This was done to prevent "dark 
side" effects (such as saddle points) , from masking hot spot 
effects. The "hot spot" consisted of biasing the potentials 
or. nine surface cells to 5000 volts negative with respect 
to spacecraft ground. Since NASCAP linearly averages surface 

cell potentials to get nodal potentials, the effective area 

2 2 

of the hot spot was 4 AX or 0.04 m . The radius of a disk 
with the same area would be 0.113 meters. Figure 9.2 shows 
the test object. In Table 9.2 we show the surface cell 
potentials along line x = 0 when V Q had floated to -5880 and 
compared to 


V 


5880 


theory 


arcsm 


(oop) 


(9.3) 


Not surprisinaly , the comparison is quite good, much 
better than 10 percent. This case was run primarily as a 
test of usefulness of the particular NASCAP object for repre- 
senting a circular hot spot, as well as to demonstrate the 
validity of the simple insulating theory. 
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TABLE 9.2. PERFECTLY INSULATING HOT SPOT POTENTIALS 



V » 5880 
o 

O’ 

It 

o 

• 

H* 

R 

v 

theory 

^N AS CAP 

.2 

-2248 

-2110 

.3 

-1446 

-1480 

.4 

-1072 

-1070 

.5 

-843 

-838 

.6 

-709 

-688 


9.2 PHOTOCONDUCTIVITY EFFECTS 

The low energy electrons reflected by the potential 
barrier in front of the plate can carry current in the plane. 

The effect of this is to discharge the plate somewhat from 
the perfectly insulating case discussed in the previous sec- 
tion. However, from our experience from dark side saddle- 
points we do not expect a drastic modification of the results. 

The fundamental assumption is that the potential change, no 
matter how small, is large enough to prevent all low energy 
electrons from leaving the surface. As a result all of the 
incident plasma current must be transported away from the hot 
spot region as skin current to outside of the disk. 

As a model problem we examine the same disk charged 
to potential V Q on an insulated flat plate. The disk has a 
radius a. The boundary condition on the plate is V(°°) = 0 
and the skin current, K, is zero at r = a. Current continuity 
on the plate is the fundamental equation for this system 

V • K + J =0 (9.4) 

o 

where J is the incident plasma current. Photosheath conduc- 
o 

tion is modeled by an effective conductivity 
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K - oE 


(9.5) 


Using the effective photosheath formulation derived 

f "71 

by Mandell, et al. , J we retain the inverse square depen- 
dence of the conductivity on the surface normal electric 
field, e.g., 

a « A/e£ (9.6) 

At a radius r the surface current integrated around 

the circumference must equal the integral of the incident 

current J over the disk from radius, a, of the hot spot to 
o 

r • 


K2irr = 


r 

f J o (r,) 
0 


2ur 'dr ' 


(9.7) 


For simplicity let us assume J Q constant. Since we are assum- 
ing photoelectrons do not escape, a similar assumption con- 
cerning secondaries is equally valid. This helps justify J Q 
being a constant. 


K2irr = J tt 
o 




(9.8) 


Substituting for K 


2 2 

2^roE | | = J q tt (r - a ) 


A E< 


2rrr 


^ — 

1 


t , 2 2, 

J Tr(r - a ) 
o 


(9.9) 

(9.10) 


We will expand around the zero conductivity solution, assuming 
small perturbations. Since for the zero conducting case 
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(9.11) 



0 , we get 



2Ar 

T 


J o (r - a > 


where the superscript zero identifies the zero conductivity 

To relate 


case and the one indicates the perturbed solution. 

Ej^ to a AV(r) i V^fr) - V°(r) we must estimate the effective 
capacitance per unit area of the surface as a function of the 
radius r. To do this we have made the assumption that the 
capacitance is the same as the mean capacitance per unit area 
of an isolated disk of radius r. First we calculate the 
charge per unit area of a disk at potential V. At large radii 
(d >> r) the potential looks like 


, _ 2 Vr _ Q disk 

^ Trd~ 4 7T€ d 
o 



(9.12) 

(9.13) 

(9.14) 


2 

The term 2rrr is the area of both sides of the disk, which 
enters in tb*' perfectly insulating case. 

Substituting we get a simple algebraic expression for 
AV in terms of E? > 


2 3 . _o 

*r r A E, : 


( AV) = 


8 (r 2 - 


a 2 ) J. 


(9.15) 
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but 



Substituting for E? 


(AV) 2 


V 





A 


From Reference 7 


(9.16) 


(9.17) 


A » 4 <e> v' 1 


sec 


(AV) 


ttV 


a r 2 J sec 

37? <E> TT 


2 2 
(r - a ) 


(9.18) 


To examine the usefulness of this expression we ran the iden- 
tical NASCAP case in the insulating case, but included secondary 
sheath conductivity as per Reference 7. The relevant values 
are 

<e> = 2 eV 

7 . 6 

o 

V = -5880 
o 

and the effective 
a * 0.1128 
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was used in the perfectly insulating case. The comparison of 
theory and experiment ia shown in Table 9.3. While the agree- 
ment is not as good as in the perfectly insulating case, it is 
still better than 20 percent. 

TABLE 9.3. HOT SPOT WITH PHOTOCONDUCTIVITIES 


■ -5890 a - 0.1128 

o 


R 

^ V theory 

0 

AV 

av NASCAP 

v 

theory 

V NASCAP 

.2 

531 

570 

-1717 

-1540 

.3 

250 

390 

-1196 

-1090 

.4 

169 

282 

-903 

-788 

.5 

129 

174 

-714 

-664 

.6 

105 

117 

-604 

-571 


9.3 CONCLUSIONS 

Both the simple theory and NASCAP predict hot spots 
to affect insulating surfaces on a scale of the hot spot 
size. The similarity between hot spot effects and other 
potential barriers seen in the sunlit sphere and SCATHA cal- 
culations is more than coincidental. The dominant physical 
mechanism is the photo and secondary electron limiting which 
sets small normal electric field boundary conditions on the 
Poisson problem. Thj s makes the exposed insulator seem to 
"disappear". The photosheath conductivity is effective in 
diminishing hot spot effects, but its inverse square de- 
pendence on normal electric field strength minimizes any dif- 
ference in the solution. 

It is worth noting that not only was the conducting 
sheath in the NASCAP example useful as a test of the simple 
theory, but it was also a good test of the NASCAP photo- 
sheath conduction algorithm. The sensitivity of the conduc- 
tivity to small field changes led to small ("-50 volt) oscil- 
lations in the final solution. 
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10. THE STRUCTURE OF THE LOW ENERGY PHOTOELECTRON 

SPACE CHARGE SHEATH 

One aspect of spacecraft charging that had not been 
thoroughly investigated using NASCAP has been the structure 
of the space charge barrier formed by photoelectrons. This 
barrier, while small in voltage compared to substorm induced 
charging, is of interest both theoretically and experimentally. 
Experiments, such as SC>10 and SC-2 can make inferences about 
the sheath structure. The extent and magnitude of the sheath 
can be important in unfolding ambient electric field data. 

Previous NASCAP calculations did not take the space 
charge of the photosheath into account because the voltage 
perturbations of the sheath are extremely small compared to 
the kilovolt/meter fields set up by differential or space- 
craft charging during substorms. However, for the case of 
the unfolding of SC-10 during quiescent conditions, the space 
charge sheath is the dominant source of field differentials 
between the two halves of the dipole. This experiment has 
provided electric field strength as a function of distance 
in front of SCATHA. 

10.1 CODE MODIFICATIONS 

In order to predict the space charge barrie.r, a slight 
modification of the NASCAP explicit photosheath treatment was 
necessary. The "SHEATH" option in NASCAP predicts photocharge 
densities, but does not use them for potential calculations. 

The space charge calculations run for SCATHA were made using 
the calculated space charge, and the potentials and charge 
densities were iterated on until self-consistency was ob- 
tained between the potentials and the charge densities ob- 
tained from particle tracking. This procedure converged 
quite rapidly. Another improvement in the "SHEATH" routine, 
that is, emission of photoe.lectrons at several angles rather 
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than just normal to the surface, was made to calculate more 

accurate charge densities. * 

These modified SHEATH routines are not designed for 
general use with NASCAP for several reasons. First, the self* 
consistent routines are presently implemented only for fixed 
spacecraft ground potentials, since the net photosheath cur- 
rents are not available to the LONGTIMESTEP features of the 
code. Secondly, to track particles and iterate on the poten- 
tials would make this procedure prohibitive to use on a pro- 
duction basis because of computer time requirements. However, 
for those interested in scientific investigations of the low 
energy sheath structure the keyword in the RDOPT file is 
"SHEATH SELF CONSISTENT". A number can be included on the 
card after these keywords to specify the charge density re- 
laxation parameter, a; 


ap + 

new 


( i-a) p 


old 


0 7 . 1 ) 


A value of unity for a would correspond to explicit iteration. 
The default value is 


^default 


0.5 


( 10 . 2 ) 


With the default value of the relaxation parameter a typical 
run converged to within ±5 percent in the charge density 
after five iterations. 


10.2 RESULTS 

Tests of the new SHEATH routines were performed using 
a simplified version of the SCATHA one-grid model. A zone 
size of 0.23 m was used on a SCATHA model with no booms. The 
surfaces were all treated as conductors with a photoyield of 
2 ‘ 10 ^ A/m^. Calculations were performed at fixed satel- 
lite potentials of +0.5, +1.0, and +5.0 volts. Figures 10.1 
tnrough 10.6 illustrate the results for the most interesting 
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Figure 10.1. Self-consistent sheath contours around a simpli- 
fied SCATHA model , vehicle potential =+0.5 volts 
top view. Contours are in units of code units of 
charge per cubic mesh unit, from 0. to -0.24 in 
steps of 0.02. The zone size is 0.23 m; to con- 
vert the contour levels to coul/nH, multiply by 
e o /(0.23) 2 = 1.67 x 10" 10 . 







7 -t -3 -1 i 3 C 7 9 It 13 IS 17 19 II S3 tS 

Y-AWf 


Figure 10.3. Photoelectron trajectories from bellyband cells 
for simplified SCATHA model, vehicle potential = 
t0. 5 volts. Lower energy electrons are returned 
to the vehicle by the space charge barrier. 
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Figure 10.4. Photoelectron trajectories from bellyband cells 
for simplified SCATHA model, vehicle potential = 
+0.5 volts. Lower energy electrons are returned 
to the vehicle by the space charge barrier. 
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Figure 10.5. Self-consistent potential contours around a sim- 
plified SCATHA model, vehicle potential * +0.5 
volts, top view. Contour spacing * 0.05 volts. 
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I Figure 10.6. Self-consistent potential contours around a sim- 

j plified SCATHA model, vehicle potential * +0.5 

volts, side view. Contour spacing = 0.05 volts. 
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case, when the vehicle is at +0.5 volts. The space charge 
density is shown in Figures 10.1 and 10.2. The irregular 
appearance of the contours is a reflection of the rather 
small number of particles tracked, five energies at each of 
five angles for every surface cell. Typical particle tra- 
jectories are shown in Figures 10.3 and 10.4; the effect of 
the space charge barrier is clear for the lower energy parti- 
cles. A barrier of approximately 0.75 volts forms at a 
distance of 75 cm above the emitting surfaces, as shown in 
Figures 10.5 and 10.6, leading to a potential minimum at 
-0.25 volts. 

Similar calculations at +1.0 and +5.0 volts satellite 
potential are shown in Figures 10.7 through 10.10. The space 
charge perturbation is smaller in these cases, so that in 
the +5.0 volt case the space potentials are everywhere posi- 
tive, the only effect being a slight compression of the con- 
tours on the sunlit side. 
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Figure 10.7. Self-consistent sheath contours around a simpli- 
fied SCATHA model, vehicle potential * +1.0 volt, 
top view. Contour levels are the same as in Fig- 
ure 10.1. 
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Figure 10.8. Self-consistent potential contours around a sim 
plified SCATHA model, vehicle potential * +1.0 
volt, top view. Contour spacing >* 0.1 volt. 
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Figure 10.9. Self-consistent sheath contours around a simpli- 
fied SCATHA model, vehicle potential « +3 volts, 
top view. Contour levels are the same as in 
Figure 10.1. Multiple maxima reflect the use of 
discrete energies to model a Maxwellian energy 
distribution. 
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Figure 10.10. Self-con3istent potential contours around a sim- 
plified SCATHA model, vehicle potential * +5.0 
volts, top view. Contour spacing * 0.5 volt. 
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11. SCATHA BODY SHADOWING 


Shadowing tables for individual SCATHA experiments 
have been described previously. J For considerations of 
the overall response of the spacecraft to the environment, 
similar tables for the spacecraft body will be useful. This 
chapter describes a set of tables and formulae which have 
been developed to provide shadowing information for body 
elements of the SCATHA satellite. Section 1.1.1 describes the 
general features of the body shadowing problem, and Section 
11.2 describes the set of tables and formulae which have been 
developed. 

11.1 BODY SHADOWING PROBLEM DEFINITION 

The correspondence between the rectangular coordinate 
system of the SCATHA satellite and the spherical coordinate 
system used to define the sun direction vector is illustrated 
in the diagram below. 

The Y-Z plane is defined by the delta attach fitting 
interface (about 5.1 inches below the bottom of the substrate). 
The sell boom lies along the negative Y axis ($> ■ 180°, 

9 - go 0 ) . The top of the satellite (where the SC1-3 is lo- 
cated) is at X » 73.82 inches while the bottom is at X » 5.10 
inches. In orbital configuration 85° ^ 0 <_ 90®. (The same 
coordinate systems were employed in Reference 8.) 

Imagine dividing the exposed SCATHA body into distinct 
elements, as shown in Figure 11.1. 

• forward surface 

• aft surface 

• top solar array 

• bottom solar array 

• bellyband 
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S ■ direction vector to sun 

9 ■ polar angle (angle between +X axis and S) 

♦ a azimuthal angle (angle between +Y axis and pro- 
jection of S onto Y-Z plane) 

+X axis is the cylinder axis of the main satellite body 


The problem is to calculate effective illuminated areas, A^, 
for each element as a function of solar angles 9 and $ . The 
total solar energy incident upon the i fc ^ element, E^, is then 
given by 


E i (9 ,$) 


A.(9,<J>) • I 


where I is the solar intensity (energy/area) . Note that since 
the spacecraft body is a cylinder of radius 33.6 inches and 
height 68.7 inches, the total solar energy incident on the 



body is approximately E * 4617 I where I is in energy/sq. in. 
Because of their relatively large lengths, the boom surfaces 
will collect significant solar energy, although the importance 
of such collection is minimized due to the weak capacitive 
coupling of these surfaces to spacecraft ground. 

11.2 SHADOWING TABLES FOR BODY ELEMENTS 

The shadowing information required for body elements is 
to ba used in analyses of the overall response of the vehicle 
to the spacecraft environment. The accuracy required for such 
considerations is much less than that needed for analyses of 
individual experiments. Accordingly, reasonable approximations 
have been made to simplify the calculations where appropriate 
for the various body elements, as discussed below. None of 
the approximations employed will lead to errors of more than 
5 percent in calculated effective illuminated areas. 

B ellyband Surface 

Shadowing of the bellyband surface involved the most 
effort since the short booms cast shadows mainly in this region. 
The five short booms were treated as 1.7 inch radius cylinders. 
The 3.5 inches radius spheres at the end of the SC2 booms were 
also distinguished as shadowing elements. Experiments mounted 
at the ends of the remaining booms were not distinguished, but 
the lengths of these booms were adjusted slightly to approxi- 
mate the shadowing effects. The cylindrical shadowing ele- 
ments used are listed in Table 11.1. With these approximations, 
the projections of the bellyband and shadowing elements can be 
analytically projected into two dimensions, and the effective 
illuminated areas calculated directly. The results are dis- 
played in Table 11.2, in 1° increments for polar angle 9 and 5° 
increments for azimuthal angle 4> . The effective illuminated 
area of the bellyband averaged over a full rotation is given 
in Table 11.2 for values of the polar angle 9 . 
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TABLE 11.1. CYLINDRICAL SHADOWING ELEMENTS FOR BELLYBAND 


Boom 


Length 


Azimuthal Angle 


SCI 97 inches 343° 
SC2-1 118 inches 303° 
SC2-2 118 inches 123° 
SC6 124 inches 12° 
SC11 164 inches 180° 


The SCI# SC6, and SC11 boom elements are slightly 
longer than the booms themselves to approximate 
shadowing by the experiments at the boom tips. 

All cylinders have diameters of 1.7 inches. 


Top/Bottom Solar Arrays, Aft Surface 

The top and bottom solar arrays and the aft surface 
are not significantly shadowed. The largest shadowing of 
these three areas occurs for the bottom solar array, where 
the SC2 booms can shadow up to 4 percent of the maximum ef- 
fective illuminated area. Referring to the dimensions in 
Figure 11.1, the effective areas are then as follows: 

Top solar array, A = 1835 sin9 

Bottom solar array, A = 1976 sin9 

Aft surface, A ■ 2425 sin (s - , 0 > tt/2 

A = 0 , 9 < it/2 

(all in square inches) . 

Forward Surface 

The lower portion of the OMNI antenna and the SC9 ex- 
periment shadow the top surface. The NASCAP SCATHA model 
was used to represent the shadowing of the top surface, and 
the HIDCEL shadowing features of the code were employed 
directly. The error introduced by the limited resolution of 
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vBLE 11.2. EFFECTIVE ILLUMINATED AREAS FOR BELLYBAND 
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All entries in square inches. 


TABLE 11.2. EFFECTIVE ILLUMINATED AREAS FOR BELLYBAND (CONTINUED) 
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All entries in square inches. 


TABLE 11.2. EFFECTIVE ILLUMINATED AREAS FOR BELLYBAND (CONTINUED) 
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the model is no more than a few percent. Table 9.3 gives the 
results, which were renormalized to the correct total area 
presented by the cylindrical spacecraft forward surface. En- 
tries for the regions <|> < 120 and <?> > 300 are omitted, since 
the SC9 cluster does not shadow the top surface in this 
region, and the OMNI shadowing is essentially constant. 
Therefore the missing entries are identical to the * 120 
values . 


11.3. 

EFFECTIVE 

ILLUMINATED 

AREAS FOR 

TOP SURFACE 

PHI 

85 


THETA 




86 

87 

88 

89 

120 

286. 

229. 

172. 

114. 

57. 

130 

280. 

224. 

168. 

112. 

56. 

1«»0 

273. 

218. 

164. 

109. 

55. 

150 

262. 

21C. 

157. 

105. 

52. 

160 

250. 

200. 

150. 

100. 

5U. 

1 70 

242. 

194. 

145. 

97. 

48. 

180 

236. 

191. 

143. 

95. 

48. 

190 

231. 

185. 

139. 

93. 

46. 

200 

232. 

185. 

139, 

93. 

4 6 • 

210 

234. 

186. 

1 4 Q« 

93. 

47. 

220 

236. 

190. 

142. 

95. 

47. 

230 

232. 

164. 

138. 

92. 

46. 

240 

240. 

192. 

144. 

96. 

46. 

250 

248. 

199. 

149. 

100. 

50. 

260 

260. 

208. 

156. 

104. 

52. 

270 

275. 

22 0. 

165. 

110. 

55. 

260 

282. 

226. 

169. 

113. 

56. 

290 

287. 

229. 

172. 

115. 

57. 

300 

289. 

231. 

173. 

116. 

58. 

AVG 

270. 

216. 

162. 

108. 

54. 


All entries in square inches. Averages include the values 
for 4> > 300 and <1> < 120, which are identical to the <}> = 120 
values . 
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12. CHARGING OF LARGE SPACE STRUCTURES 


With the advent of the shuttle era, there is in- 
creased interest in the charging characteristics of larger 
space structures, particularly in polar earth orbit. We 
have investigator, the charging of a large sphere subject to 
the environment )..y the shuttle orbiter as it passes through 
the auroral regions in its low polar earth orbit. 

The shuttle orbiter, passing through the ionosphere at 
altitudes of a few hundred kilometers, develops electrical 
potentials through accretion of charge from the natural en- 
vironment. Under normal ambient conditions the particle 
energies viewed from the satellite range from a few tenths 
of an electron volt to a few volts. Thus, the magnitude of 
vehicle potentials are at most a few volts. However, while 
passing through polar latitudes the vehicle may be subjected 
to a substantial flux of energetic electrons moving through 
the auroral zone following their injection in the magneto- 
sphere. This may cause charging to high potentials. 

Most experimental studies of spacecraft charging in 
low earth orbit have concerned small objects ('x-l m) moving 
through the ionosphere. In the absence of energetic precipi- 
tating electrons, the magnitude of the observed electric 
potentials on the INJUN 5 satellite were less than a few 
volts, in accordance with theoretical expectations.^ 17 ^ 

Even during impulsive precipitation events, observed poten- 
tials did not exceed -40 volts negative. 

More recently, theoretical studies have focused on 

charging of large objects. Parker has presented a method 

for computing sheath structures of large spherical bodies 

with high-voltage surfaces and with photoelectric/secondary 
f 181 

emission. McCoy et al . have considered problems associ- 

ated with the operation of large, high-voltage solar arrays 
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f 191 

in the ionosphere. 1 Liemohn has considered the electrical 

charging of the shuttle orbiter in the absence of fluxes of 

energetic precipitating electrons. Inouye et al. 

investigated the charging of a space based radar system 

f 211 

having an antenna with a diameter of about 70 meters. * 

Their calculation of electrical potentials in the presence 

of energetic particles are based on the application of orbit 

limited theory of Langmuir and Mott-Smith to determine the 

r 22 1 

currents of attracted species . 1 J 

The investigations of charging presented below are for 
the regime where body dimensions are large compared to the 
relevant Debye length. In this regime the currents of at- 
tracted species are estimated by adapting the large spherical 

probe theories of Langmuir and Blodgett^ 23 * and Al’pert 
f 24 1 

et al. 1 We examine the charging of a conducting sphere 
subjected to intense fluxes of energetic electrons. Factors 
relevant to a more thorough analysis of complex objects with 
dielectric surfaces are summarized. Conclusions are given in 
the final section of this chapter. 

12.1 ANALYSIS 

The purpose of the following analysis is to estimate 

the magnitudes of potential that develop on objects in low 

earth orbit (200 to 400 km) when subjected to high fluxes 
2 

(^200 uA/m ) of hot (5 to 10 keV) precipitating magneto- 
spheric electrons. Nominal values of the satellite and en- 
vironmental parameters relevant to the analysis are sum- 
marized in Table 12.1. 

We are concerned primarily with the possibly large 
negative potentials that may oe produced by the currents of 
hot electrons incident from the magnetosphere. Questions 
related to the satellite wake and its structure are not con- 

• 3 

sidered; we consider the ram ion current density NeV ^ 10 

2 ° 
amp/cm apparent to a co-moving observer as the only relevant 
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TABLE 12.1. NOMINAL VALUES OP PARAMETEFS WHICH INFLUENCE 
ELECTRICAL CHARGING IN LOW EARTH ORBIT 


Sphere Diamatar 

1000 cm 

Satallita Valocity V Q 

8 x 10 5 cm/sac 

Ambient loa Tamparatura 9 ^ 

.1 - .5 aV 

Ambiant Elaetron Tamparatura 9 # 

.1 - .5 eV 

Precipitating (Hot) Elaetron 
Tamparatura 8 p 

5-10 kaV 

Nautral Atom Density (0) 

10 10 ea* 3 

Ion Danaity (0*) 

10* - 10® cm"* 

Ambiant Oabya Langth 

<1 cm 

Thermal Electron Laraor Radius 

2 cm 

Hot Electron Larmor Radius 

400 cm 

Ion Larmor Radius 

300 cm 

2 

Currant Density (amp/cm ) 


Thermal Electron j 

10" 7 

Thermal Ion (0 + ) 

10" 10 

Photoalectron j 

io” 9 

Pracipitating (Hot) 
Electron j. 

P 

2 x 10" 8 

F 

Ram Ion j r 

CO 

1 

o 
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attribute of the satellite motion. Thus, for example, it 
is anticipated that the V Q xB inductive electric fields are 
small relative to the electrostatic fields produced by 
charging. 

To proceed further, let us first neglect the magnetic 
field. The effects of a magnetic field Vi.ll be discussed 
later. The flux of hot electrons to the satellite is assumed 
unidirectional. Since the ram ion energy (E Q * 5 eV) is 
much larger than the ion temperature, the ram ion flux will 
also be considered unidirectional. In the absence of elec- 
tric potential the precipitating electron and ram ion currents 

2 2 

to the satellite will be j irR^ and j*R . respectively. 

p o r o 

For negative potentials electrons are repelled and the 
current of precipitating electrons at the satellite is ap- 
proximately j irR exp{e<j>/0 }. This is an accurate approxi- 

r u r 

mation if the effective collection radius R is not much 

c 

greater than R . that is, if the thickness R - R of the 

O CO 

non-neutral space charge region around the object is less 
than the satellite radius. For all practical purposes in 
the cases of interest, -e<J> >> © e . Thus the cold plasma 
electrons do not enter the sheath regie n 

The effect of space charge upon current collection in 
low earth orbit by large high voltage objects is well-known, 
having been studied both theoretically and with laboratory 
experiments. Space charge effects dramatically reduce the 
current collected per unit area compared to those predicted 
by orbit limited theory. The I-V characteristics of a 
spherical probe with a ratio of radius to Debye length of 
10 is shown in Figure 12.1. The current collected per unit 
area at large voltages is substantially less than the very 
large Debye length orbit limited theory would predict. 

However, the auroral electron fluxes in polar earth orbit 
are incident currents which may be substantially larger than 
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^ R 


POTENTIAL (e$/8) 


Figure 12.1. 


The I-v characteristic for a spherical probe 
in a small Debye length plasma. Note how even 
at large potentials the probe collects just a 
few times the plasma thermal current. The 
dashed line is for long Debye length orbit 
limited collection. It is not applicable to 
large objects in low earth orbit. 




the ram ion currents. We are then interested in the inverse 
function , that is, the V-I characteristic (Figure 12.2). 

Note how dramatically the probe voltage must rise to increase 
the current collected per unit area. It is this steep V-I 
characteristic which forms the basis of the following analysis. 

The theory of the sheath surrounding a large spherical 

probe with radius R >> X_ at high potential |e$| >> 0 , 

o u e f231 

9 , in an isotropic plasma is given in Langmuir and Blodgett 1 1 
* r 24 1 

and Al'pert et al. 1 1 The effective collection radius R c 

for the case of ion attraction can be expressed as 

where 6 is the temperature of the attracted species and X the 
Debye length. F is an increasing function of its argument and 
hence of the satellite potential. 

In order to adapt the Langmuir-Blodgett theory as an 
approximation to the case of streaming ions, we relate the 
temperature e to the kinetic energy E Q of ions relative to 
the satellite by requiring that current entering the sheath 
in the isotropic and streaming cases be the same, 

N V R^ » 4itA 2 N (80AM) 1/2 (12.2) 


giving 


uMV* 

• --r 5 

where M is the 
X * 74 3 



ionic mass. 
(N/0) 1/2 cm 


The equivalent Debye length is 


(12.3) 


(12.4) 


204 


POTENTIAL (e*/0) 

10.0 30.0 30.0 40.0 90.0 


v-i 



0.0 10.0 30.0 30.0 

CURRENT <j/j 




V X D - 0 

(ORBIT LIMITED) 


T” 

40.0 


th 
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Figure 12.2. The V-I characteristic for a spherical probe in 
a small Debye length plasma. Note how even a 
small increase in probe current causes a very 
large change in the potential of the sphere. 

The dashed line is for long Debye length, 
orbit limited collection. 


Table 12.2 gives values of R /R as a function of z = (e$/0) 
4/3 ® ° 

(X/R ) ' . For values of R /R^ < 1.05, the collection 

O CO — 

radius and potential are related by the plane electrode 
Child-Langmuir law 



1 + 


2/7 ,3/4 
T~ Z 


(12.5) 


with an accuracy better than 3 percent. 

The potential on the sphere is determined by balance 
of currents. 


2 e * /6 p 2 

* R o 3 p (1 ’ s p )a = * R c j r {1+s i ) + X v 


( 12 . 6 ) 


where Sp(s^) is the total secondary yield from electron (ion) 
impact and 1^ is the total photoemission current. 

Defining 


t p - yi-v 


j = j (1+s.) 
J r J r i 


as effective electron and ion current densities corrected for 
secondary emission, Eq. (12.6) becomes 

/ R \ 2 I 

Figure 12.3 shows the dark potential on spheres of 0.5 
and 5 m radius as a function of ratio of precipitating elec- 
tron to ram ion current densities in a plasma with ambient 
5 -3 

density 10 cm . For a given current ratio the potential 
on the sphere scales roughly as the radius. More precisely, 
the potential scales with radius as (R 0 />)^ 2 for |e<J>| << 0p, 
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TABLE 12.2. EFFECTIVE COLLECTION RADIUS AS FUNCTION OF 
z = ed>/0(X/R o ) 4/3 


V R P 

1.00S 

1.018 

1.030 

1.050 

1.100 

1.150 

1.200 

1.250 

1.300 
1.340 

1.400 
1.450 

1.500 
1.600 

1.700 
1.800 

1.900 
2.000 
2.100 
2.200 

2.300 

2.400 

2.500 
2.600 

2.700 
2.800 

2.900 


z 

.001 

.005 

.010 

.019 

.052 

.094 

.143 

.199 

.264 

.337 

.421 

.510 

.610 

.833 

1.092 

1.384 

1.711 

2.074 

2.479 

2.919 
3.400 

3.920 

4.479 
5.113 
5.752 
6.472 
7.196 


Satellite Potential (volts) 



Current Density Ratio 'Tp/j* r ) 

Figure 12.3. Satellite potential as function of current 
density ratio. 
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but somewhat more slowly with R Q /X as | e<*> | increases. Ob- 
serve that the potential is an extremely sensitive function 

of J /T for values of this ratio near unity, especially for 
P ^ 

the larger sphere. 

The theory predicts that the 5 m sphere will charge 
to about the 1 kilovolt level for electron to ram ion current 
density ratios of only about two. This is to be contrasted 
with the result predicted by orbit limited ion collection. 

The approximate dark current balance 


h . (i ♦ 


e$ /e. 


e / 


for orbit limited collection predicts, for example, that 
j /j r = 300 would be required to sustain a 1 kilovolt poten- 
tial on the sphere. 


12.2 DISCUSSION 


Several effects have been neglected in determining 
that hot electrons precipitating from the magnetosphere can 
charge a large object to kilovolt potentials. We shall now 
argue that accounting for these effects will not alter the 
conclusion that such high potentials should be expected for 
the assumed charging environment. 

Consider first the effect of a magnetic field on the 
ram ions entering the sheath surrounding the satellite. A 
component of magnetic field perpendicular to the satellite 
velocity will tend to insulate the surface from the ram ion 
currents, leading to larger negative potential of the satel- 
lite. For cases of interest however, the effect is negli- 
gible. A measure of the size of this effect is given by 


a 


1 M ,s 2 A 2 

? M “ ci d 2 2 

* Z w . T 

ci 

where M is the ion mass, w ci its gyrofrequeney, d the thick- 
ness of the sheath, and t the flight time of an ion across 
the sheath. For the cases represented in Figure 12.1, d < R, 
so that 


a < 0.2/| e<j> (volts) | 

which is negligibly small except at very low levels of satel- 
lite potential. 

The hot electrons responsible for charging the satel- 
lite were considered to approach the space charge sheath uni- 
directionally , as pertains in the limit of strong magnetic 
fields where the Larmor radius is small compared with the 
radius of the satellite. More probably, the electrons, be- 
cause of their pitch angle distribution, would enter the 
repulsive sheath with a more nearly isotropic distribution of 
directions. Assuming that the one sided thermal plasma cur- 
rent densities are the same in the unidirectional and iso- 
tropic limits, the effective electron current toward one hemi- 
sphere of the satellite in the isotropic limit is twice that 
which pertains in the unidirectional case. In the absence 
of no other effect associated with the magnetic field, the 
result would be greater charging. 

The charging current given by Eq. (12.6) for the case 
of repelled electrons incident unidirectionally from infinity 
applies in the limit of zero gyroradius. In the opposite 
limit of vanishing magnetic fields, again assuming that 
electrons enter the sheath unidirectionally, fewer electrons 
reach the satellite because of the deflection by the repul- 
sive electric field. The reduction in current is small how- 
ever, and the charging current accurately represented by 
Eq. (12.6) provided that the repulsive potential on the 
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2 

satellite satisfies (e<j>/0 p ) << 1. This requirement, which 

is satisfied in the case of Figure 1 for potentials less than 
about 2 kV, follows from the conservation laws of energy and 
angular momentum which permit one to express the current to 
the satellite as 


I 



Essentially, the electron current crossing the sheath is not 
substantially modified by the magnetic field, a circumstance 
we expect to pertain is long as 


( Uee T) 2 £ m d J /6 


<< 


1 


This condition is well satisfied for potentials in Figure 12.3 
at the kilovolt level. 

In applying the Langmuir-Blodgett probe theory, we have 
neglected the contribution of electrons to the space charge in 
the sheath. This is a valid approximation because the velocity 
of electrons in the sheath is large compared with ion veloc- 
ities, except perhaps for the contribution of secondary and 
photoelectrons near the surface. Near the surface, however, 
the electric fields are dominated by surface charge and little 
affected by space charge. 

Secondary and photoelectrons move through the sheath 
with smaller energies than the precipitating magnetospheric 
electrons and are therefore more strongly affected by the 
magnetic field. The potential developed by the satellite is 
affected however only if the emitted electrons return to the 
surface, leading to higher potentials than if the electrons 
escape. 
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In all previous considerations, we have supposed that 
the satellite is a conducting sphere. The shuttle orbiter is 
actually a geometrically complex object whose surface is 
coated with dielectric materials, and both ion and electron 
fluxes are apt to be strongly heterogeneous functions over 
the satellite's surface. The degree of heterogeneity will 
be affected by the geometry of the satellite, its motion 
through the ionosphere, the variation of surface properties, 
such as secondary yield, and by the magnetic field. Un- 
doubtedly the sheath surrounding the orbiter will have a 
complicated geometrical structure not easily represented by 
simple spherical probe models. Multidimensional computer 
models will be required to determine the strong differential 
voltages which are expected to develop on the vehicle. 

12.3 CONCLUSIONS 

Ambient currents of hot electrons (5-10 keV) of 
2 

200 uA/m will charge a 5 meter sphere in low polar earth 
orbit to kilovolt potentials in eclipse. Such potentials 
are about 1 order of magnitude larger than occur for smaller 
satellites (^R c % 0.5 m) in a similar orbit. On this basis, 
one should expect negative potentials of around 1 kilovolt 
to develop on the shuttle orbiter. Because of the dielectric 
coating on the orbiter, and the non-uniform character of the 
charged particle fluxes expected at the vehicle's surface, 
differential surface potentials of the order of one kilovolt 
should also occur. 
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13. CONCLUSIONS AND RECOMMENDATIONS 


The simulations of the charging response of the 
SCATHA spacecraft described in the above chapters were based 
on the actual observed plasma environments as well as proto- 
typical representations of the geosynchronous ambient plasma. 
This chapter provides a brief summary of the major simulation 
results, and gives an assessment of the validity and implica- 
tions of these results when possible. The studies to date, 
together with preliminary reports of data from the SCATHA 
spacecraft, suggest several avenues for extension of the 
simulations, and these avenues are briefly described as well. 

Perhaps the most important result of the simulations 
was a prediction originally made before the launch of the 
vehicle, ^ and emphasized by the simulations described in 
Chapter 5: the SCATHA booms were isolated as the most 

likely sites for discharges to occur. In fact, the onxy 
charging related malfunction observed, the SC2 experiment 
failure and the associated telemetry upset on day 89, seems 
to have been caused by just such a discharge. Continuing 
analysis of this event is in progress, but the observed ef- 
fects are consistent with a scenario based on a surface 
f lashover along the SC2 booms. Thus the NASCAP/SCATHA model 
served one of its major purposes: likely discharge sites on 

an operational spacecraft were correctly identified in ad- 
vance of their occurrence. 

Simulation of the "passive" charging event that oc- 
curred during eclipse on day 87, 1979, described in Chapter 6, 
represents the first unambiguous test of the NASCAP code. 

It is a test that was passed most convincingly. These re- 
sults, coupled with the realistic potentials predicted by 
the test environment simulations reported in Chapter 5, 
allows us to conclude that NASCAP is a predictive code, and 
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that we can calculate the potentials of spacecraft under 
realistic conditions with added confidence. 

Simulations of charging in sunlight have revealed two 
interesting features. The first effect, which was also ob- 
served during the eclipse charging in the moderate temperature 
plasma, was the importance of field limiting of low energy 
fluxes (secondaries and photoelectrons) in decreasing pre- 
dicted equilibrium potentials. Such field limiting is the 
likely mechanism which led to the observed sunlight charging 

of SCATHA to a few hundred volts negative on several occa- 
r 9 1 

sions. 1 This same mechanism can also prevent the SSPM 
reference band from tracking plasma ground in sunlight. The 
second effect of interest was the persistence of differential 
potentials in the simulation of the exit from eclipse, which 
led to some spacecraft surfaces being driven to positive 
potentials for brief periods. Since the vehicle upset on 
day 89 was associated with high positive vehicle potentials, 
the transient response during eclipse exit bears further in- 
vestigation. 

Three specific design suggestions for passive mitiga- 
tion of spacecraft charging effects emerge from the studies 
described above : 

1. Eliminate exposed surfaces weakly capacitively 
coupled to spacecraft ground, as typically repre- 
sented by the SCATHA booms. 

2 . Include material charging response as a design con- 
sideration. Use of high atomic number materials 
having high backscattering yields and coatings having 
very high secondary yields (MgF2 on solar cell covers) 
along with materials having more typical charging 
characteristics guarantees the development of dif- 
ferential charging. 

3. Avoid large insulating areas which are shadowed for 
long periods, such as the white paint coating on the 
SCATHA aft surface, in order to reduce field limiting 
of low energy emission. 
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Perhaps the most serious question concerning space- 
craft charging that still remains unanswered is the mechanism 
of natural discharge. Even though the potential for damage 
arising from such discharges is a major motivation behind the 
effort to understand charging* a successful model explaining 
how and when* discharges occur has yet to be developed. In- 
cluded among the many suggestions that have been made is 
the correlation between discharge and the radiation- induced 
bulk conductivity, pointed out in Chapter 4. Future analysis 
of SCATHA charging should focus on those days when discharges 
were either unusually frequent or unexpectedly absent, and 
include the entire plasma spectrum up to MeV energies. 
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APPENDIX A 


SCATHA MODEL MATERIAL PLOTS GENERATED 
BY NASCAP CODE 





Surface cell material composition as viewed from the positive Z direction 
for Z values between 1 and 33. 







ition as viewed from the negative Z direction 
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APPENDIX B 


ANALYTICAL SOLUTION OF 1-D POISSON EQUATION 
FOR SC4 ION GUN MODEL 

For the planar case, the Poisscn equation for space 
charge limited emission of cold ions into an ambient plasma 
whose shielding effect is assumed to be proportional to the 
potential is given by (esu) 



-4tt1 


m m 


*^2 

1/2 xi 


(B.l) 


where <J> A is the potential of the virtual anode at x ■ 0 and 
j is the current density. We seek solutions of Eq. (B.l) 
which satisfy the boundary condition v 3$/dx) « 0 at x » 0 
and decrease monotonically as x increases for zero to ». 

The first integral of Eq. (B.l) gives 



8tt 



(< P a - 4>) 1/2 



G(<J>,<|> a ) 


The electric field vanishes for given by 
G( ^o' d> A ) * 0 


(B. 2) 


and the position of this zero occurs at x * « if 



(B. 3) 


This last condition is just equivalent to the vanishing of 
the right hand side of (B.l) for <J> * 4> 0 * 
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Equations (B.2) and (B.3) give 


*o 

a D 


,- 1/2 


w{lr) <V*o> 


- 1/2 


2 2 

* - <j>7 

y o y A 


It follows that 


/ 2 ,r 1/2 

‘Mir) <V*o> 


1/2 



*o * V 3 


and 


l_/2e\ 1/2 ( W 


3/2 


8 it \m / 


15 


i 


224 


c 


APPENDIX C 

FORMAT FOR TABULATED SPECTRAL DATA 

For NASCAP operating in the DIRECT mode to be able to 
read on the tabulated data, it must be prepared according to 
the following specifications and format. 

1. Magnetic Tape Characteristics 

Spectra can be provided on a coded 9-track magnetic 
tape with the following characteristics: 

unlabelled 
1600 bpi 

EBCDIC or ASCII coded 

fixed length record: 80 characters per record 

fixed block size: 20 records per block 

2 . Data Format 

Each data tape will consist of header records followed 
by repeated series of data records. The data will be read by 
a FORTRAN program using the FORMAT statements indicated below. 

HEADER RECORD 1. DETECTOR FORMAT (80A1) 

Identifies the detector (s) used to obtain the data. 

HEADER RECORD 2. SOURCE FORMAT (80A1) 

Identifies the individual (s) responsible for preparing 
the data. 

HEADER RECORDS 3 through 10. COMMENTS FORMAT (80A1) 

Any relevant information regarding the data can be 
included here, such as date data tape was generated, 
detector mode of operation, and what corrections 
have been applied to the raw data. 


f 
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HEADER RECORD 11. YEARq, DAY Q , SEC Q , 

YEAR 1# DAY 1# SECj^ FORMAT (10F8.0) 

YEARq, DAYq , SECq * time of earliest spectrum on tape 
YEAR^, DAY^, SEC^ « time of latest spectrum on tape 
Each series of data records will represent a complete 
energy scan by the detector. 

DATA RECORD 1. YEAR, DAY, SEC, NBINS, DELTA, VSAT, 

|3|, SX, SY, SZ FORMAT (10F8.0) 

YEAR, DAY, SEC * time energy scan was begun 
NBINS ■ number of distinct energy bins in the scan 
DELTA ■ time (seconds) between each data point in the 
scan of the spectrum 

VSAT * satellite potential during scan (volts) 

|$| » sun intensity (1.0 * full sun) 

SX, SY, SZ « normalized sun direction vector components 
at start of scan 

DATA RECORDS 2 through (NBINS+1) . ENERGY, log 1Q (F i ), 

log 10 (F e ) ' a ' BX ' BY ' BZ FORMAT (10F8.0) 
Each of these records represents a data point on the 
scan of the energy range 
ENERGY » energy (eV) 

F. = ion distribution function (sec 3 /m 6 ) 

1 3 6 

F g = electron distribution function (sec /m ) 

ft * detector view angle (degrees) 

a = pitch angle (degrees) 

BX, BY, BZ = magnetic field vector components 

(nT * 10~ 9 W/m 2 ) 

(The a value is redundant since it can be calculated 
from ft and the magnetic field vector.) 

DATA RECORD NBINS+2 . (END OF DATA MARKER) FORMAT (10F8.0) 
This record will contain any negative reel number to 
indicate the end of the spectral scan. (This record 
is redundant since NBINS is known.) 
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The data records 1 through (NBINS+2) are then repeated 
for each spectrum. Some of the information above may not be 
available for each spectrum or each data point. The follow- 
ing conventions can be used to indicate that the data is to 
be ignored: 

VSAT: any value greater than + 10000 

| 3 | : any negative value 

SX, SY, SZ: blank or zero 

F i' F e ! blan,c or ** to 

ft , a: any value greater than 360 

BX, BY , BZ: blank or zero 
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APPENDIX D 

FITTING THE TABULATED DATA TO FUNCTIONAL FORMS 


1. Single Maxwellian Fits 

The data was fit to a form 


f(E) • N 



,-E/T 


where m is the particle mass and T is the temperature of the 
Maxwellian. The density N is given by the zeroth moment MO 



(D.l) 


Equation (A.l) applies when the spacecraft is not charged. 
For a potential of <fi on the spacecraft the expression is 
modified . 


MO 



f (E ♦ q$) dE 


where C is the energy of the lowest energy data points in- 
cluded in the fit (i.e., the cutoff) and q is the charge on 
the particle at hand. 


f(E + q$) - f (E) 



(5 = f ) 3/2 e ^ 4/T f E 1/J e' E/T dE 

J c 

( i ¥?) V2 ^< 3 / 2 , y ) 


229 


It is easy to show that 


«° - .-«/* - .- <I * /T [.rfc(y 1/2 ) ♦ -2 

whara y - C/T. 


N - MO . , 5* /T [«fc(y l/2 ) ♦ 2 (f ) ^ .-y]' 1 
Hence we can estimate the density N by measuring the 


MOs ^ E 1/2 f(E) AE 
E-C 


The second moment, M2, has the form 


M2 


2 

T 



dE - N • T 


for am uncharged spacecraft. Introducing a cutoff C 
potential <$ leads to a result similar to that for MO: 




(D.2) 
moment MO: 

(D. 3) 

and 

(D.4) 


Equations (D.2) and (D.4) form two nonlinear simultaneous 
equations for T and N. Solution by iteration leads to values 
for N and T that make up the single Maxwellian fit. 


2 . 


Double Maxwellian Fits 


The double Maxwellian fits were made by minimizing the 
relative error (least squares) . The desired function has 
the form: 


f (E) - N x 




(D . 5 ) 


An initial choice of values for ? x and T 2 were made. The 
fit was made to agree exactly with the measured data at two 
points, one from the low and the other from the high energy 
regime. This determined the values of and Nj. and had 
the effect of weighting the fit around the fixed points and 
ensuring a good compromise fit over the whole energy range. 
All possible combinations of choices for T x and Tj, between 
realistic limits, were tried and the values that gave the 
minimum error were used as the double Maxwellian fit para- 
meters . 


3 . Discussion 

In all of the fits a cutoff of 1000 eV was used for 
the repelled species (electrons); i.e., only data above 
1000 eV was included in the fits. Using data below this value 
lead to erratic and often rather unphysical values for the 
fitting parameters. For the attracted species (ions) the 
cutoff was taken as 1000 eV or the spacecraft potential, 
which ever was the greater. In the double Maxwellian fitting 
procedure the lower limit for the choice of temperature was 
forced to be one-half of the spacecraft potential for the 
repelled species. This ensured the absence of low tempera- 
ture, high density components which were not observable an the 
surface due to the spacecraft potential. 


231 


REFERENCES 


1. Katz, I., D. E. Parks, M. J. Mandell, J. M. Harvey, 

D. H. Brownell, Jr., S. S. Wang and M. Rotenberg, 

"A Three-Dimensional Dynamic Study of Electrostatic 
Charging in Materials," NASA CR-135256, August 1977. 

2. Katz, I., J. J. Cassidy, M. J. Mandell, G. W. Schnuelle, 
P. G. Steen, D. E. Parks, M. Rotenberg and J. H. 
Alexander, "Extension, Validation, and Application 

of the NASCAP Code," NASA CR-159595, January 1979. 

3. Cassidy, J. J. , "NASCAP User’s Manual - 1978," 

NASA CR-1594 17 , August 1978. 

4. Katz, I., J. J. Cassidy, M. J. Mandell, G. W. Schnuelle, 
P. G. Steen and J. C. Roche, "The Capabilities of the 
NASA Charging Analyzer Program," Proceedings of the 
Spacecraft Charging Technology Conference, 1978, NASA 
Conference Publication 2071, AFGL-TR-79-0082 . 

5. Schnuelle, G. W. , D. E. Parks, I. Katz, M. J. Mandell, 

P. G. Steen, J. J. Cassidy and A. Rubin, "Charging 
Analysis of the SCATHA Satellite," Proceedings of the 
Spacecraft Charging Technology Conference, 1978, NASA 
Conference Publication 2071, AFGL-TR-79-0082. 

6. Schnuelle, G. W. , I. Katz, M. J. Mandell and A. G. Rubin, 
"Simulation of the Charging of the SCATHA (P78-2) Satel- 
lite," presented at the AGU Fall Meeting, San Francisco, 
December 1979. 

7. Mandell, M. J. , I. Katz and G. W. Schnuelle, "Photo- 
electron Charge Density and Transport Near Differentially 
Charged Spacecraft," IEEE Trans. Nuc . Sci., NS-26, 1979. 

8. Steen, P. G. , "SCATHA Experiment Shadowing Study," 
Systems, Science and Software Topical Report, SSS-R- 
78-3658, May 1978. 

9. Whipple, E. C. , personal communication. 

10. Mizera, P. F. , "Natural-Artificial Charging: Results 

From the Satellite Surface Potential Monitor Flown on 
P78-2," presented at AIAA Aerospace Sciences Meeting, 
January 1980. 


Ashley, J. C., C. J. Tung, V. E. Anderson and R. H. 
Ritchie, (i) AFCRL-TR-75-0583 ; (ii) RADC-TR-76-220; 

(iii) RADC-TR-76-125; (iv) IEEE Trans. Nucl. Sci., 
NS-25/6, 1566, 1978. 

Kaye, S., et al. , SC8 data. 

Frederickson, A. R., NASA Conference Publication 2071, 
AFGL-TR -79-0082, 559, 1978. 

Treadaway, M., presented at SCATHA meeting. Aerospace 
Corporation, September 1980. 

Frederickson, A. R. , IEEE Nuc. Sci., NS-24, 2532, 1975. 

Nightingale, R. W. , et al . , SC3 data, Lockheed Research 
Publication, LMSC/D7?^804. 

Sagalyn, R. C. and W. J. Burke, "INJUN 5 Observations 
of Vehicle Potential Fluctuations at 2500 km," Pro- 
ceedings of the Spacecraft Charging Technology Con- 
ference, AFGL-TR-77-0051 , NASA TMX-73537, February 1977. 

Parker, L. W. , "Plasmasheath-Photosheath Theory for 
Large High-Voltage Space Structures," Space Syste ms 
and Their Interaction with Earth’s Space Environment , 
Progress In Astronautics and Aeronautics, 71, TTT, 

1980. 

McCoy, J. E. , A. Konradi and 0. K. Garriott, ibid, 
p. 523. 

Liemohn, H. B., "Electrical Charging of Shuttle Orbiter," 
Battelle, Pacific Northwest Laboratories Report BN 5A 
518, June 1976. 

Inouye, G. T., R. L. Wax, A. Rosen and N. L. Sanders, 
"Study of Space Environment Physical Processes and 
Coupling Mechanisms," AFGL-TR-79-0206 , September 1979. 

Mott-Smith, H. M. and I. Langmuir, Phy-i. Rev., 28 , 

727, 1926. 

Langmuir, I. and K. Blodgett, Phys. Rev., 22^, 347, 

1923; 24, 99, 1924. 

Al'pert, Ya. L. , A. V. Gurevich and L. P. Pitaevskii, 
"Space Physics With Artificial Satellites," Consultants 
Bureau, New York, pp . 186-210, 1965. 


233 


